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In the limit of large number of colors Nc the nucleon consisting of Nc quarks is 
heavy, and one can treat it semiclassically, like the large-Z Thomas— Fermi atom. 
The role of the semiclassical field binding the quarks in the nucleon is played by 
the pion or chiral field; its saddle-point distribution inside the nucleon is called 
the chiral soliton. The old Skyrme model for this soliton is an over-simplification. 
One can do far better by exploiting a realistic and theoretically-motivated effective 
chiral lagrangian presented in this paper. As a result one gets not only the static 
characteristics of the nucleon in a fair accordance with the experiment (such as 
masses, magnetic moments and formfactors) but also much more detailed dynamic 
characteristics like numerous parton distributions. We review the foundations of 
the Chiral Quark-Soliton Model of the nucleon as well as its recent applications to 
parton distributions, including the recently introduced 'skewed' distributions, and 
to the nucleon wave function on the light cone. 
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1 Introduction 

Application of the QCD sum rulestl to nucleons, pioneered by B.L. lofFea has 
taught us several important lessons. One is that the physics of nucleons is 
heavily dominated by effects of the Spontaneous Chiral Symmetry Breaking 
(SCSB). One sees it from the fact that all loffe's formulae for nucleon observ- 
ables, including the nucleon mass itself, is expressed through the SCSB order 
parameter {'ijj'ip). Therefore, it is hopeless to build a realistic theory of the 
nucleon without taking into due account the SCSB. 

The chiral quark-soliton model of the nucleon we are going to review in 
this paper is based on two ideas. One is the dominating role the SCSB plays 
in the dynamics of the nucleon bound state. The second is, in fact, a technical 
tool, not a matter of principle. The idea is that our world with the number 
of colors Nc = 3 is not qualitatively different from an imaginary world with 
large number of colors. If one is able to build a nucleon as a bound state 
of Nc quarks at large Nc (keeping Nc as a free algebraic parameter) and if, 
furthermore, 1/Nc corrections are under control, one can reasonably expect 
that the picture one gets putting A^c = 3 is not too far from reality. We shall 
see that actually in some cases the leading- A^c predictions are satisfied in nature 
to a 1% accuracy (!), that is far better than one might expect beforehand from 
a generic 1/A^c ~ 30% correction. 

We start from surveying the background of the model: spontaneous chi- 
ral symmetry breaking and the effective chiral lagrangian (E^L) including the 
Wess-Zumino term. We formulate a simple but realistic E^L which is com- 
patible with the phenomenology and follows from certain theoretical consider- 
ations. The properties of this E^L is studied in some detail. For that purpose 
we introduce a technique of gradient expansion of functional determinants, 
which is widely used in many other physical applications. We then build the 
semi-classical model of the nucleon, justified at large Nc, and present another 
technical tool, that of the quantization of soliton rotations. In the context of 
the nucleon problem it is necessary to describe different spin-isospin states of a 
nucleon as well as that of the A resonance. Finally, we briefly describe recent 
development in the field, namely, the calculation of various parton distribu- 
tions inside the nucleon at low virtuality. Parton distributions obtained in the 
model satisfy all general requirements: they are positive, satisfy appropriate 
sum rules, etc. In addition, they appear to be in a good accordance with the 
data (where measured) without any fitting parameters whatsoever. 



2 How do we know chiral symmetry is spontaneously broken? 

The QCD lagrangian with Nf massless flavors is known to posses a large global 
symmetry, namely a symmetry under U{Nf) x U{Nf) independent rotations of 
left- and right-handed quark fields. This symmetry is called chiral. [The word 
was coined by Lord Kelvin in 1894 to describe molecules not superimposable on 
its mirror image.] Instead of rotating separately the 2-component Weyl spinors 
corresponding to left- and right-handed components of quark fields, one can 
make independent vector and axial U{Nf) rotations of the full 4-component 
Dirac spinors - the QCD lagrangian is invariant under these transformations 
too. 

Meanwhile, axial transformations mix states with different P-parities. Were 
that symmetry exact, one would observe parity degeneracy of all states with 
otherwise the same quantum numbers. In reality the splittings between states 
with the same quantum numbers but opposite parities are huge. For example, 
the splitting between the vector p and the axial ai meson is (1260 — 770) ~ 
500 MeV; the splitting between the nucleon and its parity partner is even 
larger: (1535 - 940) ~ 600 MeV. 

The splittings are too large to be explained by the small bare or current 
quark masses which break the chiral symmetry from the beginning. Indeed, 
the current masses of light quarks are: niu — 4 MeV, m^ ~ 7 MeV, nis — 
150 MeV. The only conclusion one can draw from these numbers is that 
the chiral symmetry of the QCD lagrangian is broken down spontaneously, 
and very strongly. Consequently, one should have light (pseudo) Goldstone 
pseudoscalar hadrons - their role is played by pions which indeed are by far 
the lightest hadrons. 

The order parameter associated with chiral symmetry breaking is the so- 
called chiral or quark condensate: 

(V^V') =i-(250MeV)3 (1) 

at the scale of a few hundred MeV. It should be noted that this quantity is 
well defined only for massless quarks, otherwise it is somewhat ambiguous. By 
definition, this is the quark Green function taken at one point; in momentum 
space it is a closed quark loop. If the quark propagator has only the 'slash' 
term, the trace over the spinor indices understood in this loop gives an identical 
zero. Therefore, chiral symmetry breaking implies that a massless (or nearly 
massless) quark develops a non-zero dynamical mass (i.e. a 'non-slash' term 
in the propagator). There are no reasons for this quantity to be a constant 
independent of the momentum; moreover, we understand that it should anyhow 
vanish at large momentum. The value of the dynamical mass at small virtuality 



can be estimated as one half of the p meson mass or one third of the nucleon 
mass, that is about 

Af (0) « 350 MeV; (2) 

this quantity is directly related to chiral symmetry breaking and should emerge 
together with the condensate (|l|) in any theoretical description of the SCSB. 

3 Effective Chiral Lagrangian (E^L) 

Let us consider QCD in the chiral limit, meaning that we put the light current 
quark masses to zero. For u, d quarks this is a good approximation to reality: 
it is known that most physical observables computed in the chiral limit differ 
from their true values by about 5%, unless the quantity is for some special 
reasons divergent in the chiral limit. For example, the nucleon mass in the 
chiral limit is lighter by about 5%: this is the contribution of the cr-term to 
the nucleon mass. For the s quark this idealization has worse accuracy, about 
20%. We would like to keep the number of light flavors Nf a free parameter 
putting it equal to 3 or 2 depending on whether we incorporate strangeness or 
not. 

Once chiral symmetry is spontaneously broken, the lightest degrees of free- 
dom in QCD are N'i — 1 pseudoscalar Goldstone bosons (we shall call them 
'pious' for short even in the case of Nf = 3 where four kaons and the rj meson 
are also Goldstone bosons). Therefore, in the low-momenta region QCD is re- 
duced to a theory of massless interacting Goldstone bosons. The appropriate 
nonlinear lagrangian describing Goldstone pious is called the effective chiral 
lagrangian (E^L). It is convenient to write the pion fields in terms of a unitary 
(iV| - 1) X (iV| - 1) matrix, 

C/(a;) = expi -^, (3) 

where A"^ are eight Gell-Mann matrices in case we consider Nf — 3 flavors 
and three Pauli matrices for Nf = 2. We divide the exponent by a pion decay 
constant -F^ = 94 MeV in order to use a properly normalized pion field with a 
standard kinetic energy, see below. 

3. 1 General properties 

In the chiral limit pious are strictly massless, so the ExL needs to have at 
least two derivatives of the matrix U , or more. Another general requirement 



for building the E^L is that it should be chiral- invariant, meaning invariance 
under a global S'C/l(^/) x SUB.iNf) rotation, 

U^AUB\ U^-^BU'^A^ (4) 

with constant SU{Nf) matrices A, B. It is useful to introduce hcrmitian ma- 
trices 

L^ = zt/a^[/t = ^a^^^A^ + ... (5) 

transforming under the rotation (^ as L^ — > ALf^A"^ . 

A general form of the E^L can be presented as a series in the number of 
derivatives of the fields; it is called the derivative expansion. The leading term 
is the two-derivative one; there is only one possibility to write it, compatible 
with the invariance (0): 

C^^^ = ^TvL^L^^^Ti-d^U^d^U (6) 

ff> ^'4\2 1 fABE fCDE q A^B-^ ^C D . nC^GN 

The first term is the usual kinetic energy for pions, while the second (quartic 
term) describes the p-wave pion scattering. The Noether current associated 
with the left-hand transformation in (|6|) by an x-dependent matrix A is 

Jt = "^ Tr {lX) - -F.d.n^ + ... (7) 

Charged pions decay owing to this current which couples to the W bosons, 
hence the normalization factor F^^ in the above equations. 
In the four-derivative order there are several invariants. 



£(4) = LiTr (£^ip)Tr (L^L,) + LsTr (L^L,)Tr (L^L,) 

+ L3Tr(i^L,,L,L,) + C4Tr(a^i^)2 (8) 

where we indroduce the coefficients Li,... in accordance with notations of 
Gasser and LeutwylerQ These coefficients have been fitted to the d-waye pion 
scattering and other low-energy data, and their present-day values areu 

Liw2L2, La = (1.35±0.3)10"^ L3 = (-3.5 ± 1.1)10"^ (9) 
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The equation Li = 2L2 is a prediction of large Nc] the term with C4 is usually 
expelled from the four-derivative family as it affects the physical amplitudes 
only at higher orders; at Nf = 2 the three first terms in Eq. (g) are not inde- 
pendent. For these reasons the four-derivative terms of the E^L are sometimes 
presented as 

£(4) = (2L2 + L3)Tr(L^L^L,L,) + L2Tr(L^i,L^L,) 

= {3L2 + L^)Ti-iLlLl) + ^Tv[L,L.f (10) 

where 2L2 + L3 is compatible with zero; the last term is called the Skyrme 
term. 

In the six-derivative order there are rather many invgujiants but only limited 
information is known about them from phenomenologyJj 

3.2 Wess-Zumino term 

There is anot|lier famous term in the four-derivative order, called the Wess- 
Zumino term.B It has several interesting peculiarities. First, it cannot be writ- 
ten down in a local form as EqaJ^, ^). One of the forms has been suggested 
by Fides and one of the authorsQ as an integral over a parameter: 

2-f,/^......Tr(e--a.e-)ygQ (11) 

This expression iz nonzero only starting from Nf > 3. The leading 0{Tr^) term 
gives the low-energy amplitude of the process KK — > tttttt. For two flavors the 
Wess-Zumino term is zero, at least perturbatively, in the sense that all terms 
of the expansion in powers of the pion field are zero. However, as shown by 
Wittenjj for certain 'large' fields it is nonzero, even for U G SU{2). Also, the 
variational derivatives of Eq. (|ll|) arc nonzero. For example, one may wish 
to 'gauge' the Wess-Zumino action by replacing the derivatives in Eq. ( PT|) by 
covariant derivatives of the flavor SUL{Nf)xSUB{Nf) group. This is necessary 
if one likes to learn how pions couple to the phpton and the W, Z bosons. 
The task of gauging Eq. (^l]) is not a trivial one El given the nonlocal form of 
the Wess-Zumino action. The complete and correct result has been given by 
Dhar, Shankar and Wadia.B The gauged Wess-Zumino action is nonzero for 
any Nf including Nf = 2. It gives, in particular, the low-momenta limit of the 
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jr^z ^ ^^ I ^^^ 



processes tt^ —> 77 and 7 —^ mnr. In a more common language these processes 
are determined by the axial anomaly. Therefore, the Wess-Zumino term in the 
ExL incorporates the axial anomaly, but in a 5J7l(-^/) x SUR{Nf) invariant 
way. 

Last, let us present the Wess-Zumino term in a 5-dim form suggested by 
WittenJj The action corresponding to the lagrangian (WJh can be written as 



S 



wz 



M ^ ^oTrTT / '^ 2;ea/3^5eTr iLaLf^L^LsL, 



(12) 



In this equation it is understood that the pion field U is continued without 
singularities from the physical 4-dim space-time to a fifth dimension such that 
the physical space-time serves as a border of the 5-dim 'disk'. Actually, it can 
be shown that the integrand in Eq. (|l2|) is a full derivative, so that Eq. (12) 
does not depend on the concrete way one continues the pion field inside the 
'disk'. The above Eq. ( pi] ) can be viewed as a specific continuation procedure, 
with the parameter s playing the role of the 'disk' radius. 

If the light current quark masses are nonzero, the pion fields become mas- 
sive (they are called pseudo-Goldstone fields in this case), and there appear 
terms in the E^L without derivatives, 



Trm(C/ + t/^), Ti mUmU\ Ti m{U + U'')L^L^, 



etc. 



(13) 



where m is the matrix of current quark masses. However, in this paper we 
shall restrict ourselves to the chiral limit and neglect such terms. 

Were QCD exactly solvable, one would be able to derive all terms of the 
derivative expansion of the E^L. The dimensionfuU coefficients (like Ft^, etc.) 
should be expressible through the only dimensional constant one has in the chi- 
ral limit of QCD, namely Aqcd, this quantity appearing, via the transmutation 
of dimensions, as a renormalization-invariant combination of the ultra-violet 
cutoff n and the gauge coupling given at this cutoff. 
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So far QCD is not exactly solved, even in the large Nc limit, therefore, to 
build a realistic E^L one has to rely upon models, physical considerations and 
comparison with phenomenology. 



8 



3.3 Skyrme model 

A very simple model is named after Skyrmell3 It consists in truncating the 
infinite series of the derivative expansion at the four-derivative level, moreover, 
leaving only one term in Eq. m), the 'Skyrme term', 

^Skyrme ^ r^^r L,L, + ^Tr [L.L^r- (15) 

Skyrme has suggested this simplified version of the chiral lagrangian to find a 
static configuration of the pion field minimizing the appropriate energy func- 
tional. He conjectured that this static configuration should be associated with 
the nucleon. More than two decades later Wittena has revitalized the original 
Skyrme's idea showing that, if supplemented by the Wess-Zumino term (13), 
a static local minimum of the E^L can be, indeed, associated with a nucleon in 
the large- A^c limit of QCD. The model of a nucleon based on thn^local minimum 
of Eq. ( p^ ) has been presented by Adkins, Nappi and WittenEJ and called the 
skyrmion. 

Simple as it is, the model cannot be fully realistic because we know that 
there are other terms with four derivatives, numerically of the same magni- 
tude as the commutator term in (|l5|), not to mention higher-derivative terms. 
Adding the Wess-Zumino term simply 'by hands' looks somewhat ad hoc, as 
is the choice of the simplified E^L. 

3.4 Low-momenta limit of QCD 

It has been first noticed in ref. Q that the Wess-Zumino term can be obtained 
as a functional integral over the quark fields. We shall show it later in section 
4 and now let us consider the interaction of quarks with the chiral field. As 
explained in section 2, an inevitable consequence of SCSB is the appearance of 
a dynamical or constituent quark mass M, a 'non-slash' term in its propagator. 
A free lagrangian, 

/;f™ = V^(i^-Af)V^, (16) 

is, however, not invariant inder axial rotations, 

t/i ^ exp(za^A^75)7/', -0 ^ -0exp(ia'^A^75), (17) 

because of the mass term. To construct a lagrangian which is invariant un- 
der axial rotations Goldstone pion fields need to be involved. The simplest 
lagrangian invariant under (O) is 

£ = ip{i^-MU^^)^, (18) 
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U^- = exp I — ^U — \-W — - — , (19) 

r.^ z z 

since the rotation of the quark fields here can be compensated by renaming of 
pion fields, U -^ AUA where A = exp(za'^A'^). 

The dynamical quark mass M is not a constant: we know that at large 
space-like quark virtualities M{k) has to vanish because of asymptotic freedom, 
to mention one reason. Therefore, the interaction of quarks with chiral fields is 
necessarily momentum-dependent, i.e. non-local. To avoid difficult questions 
of what is going on at time-like virtualities we prefer to formulate the theory 
first in the space-like region, i.e. in Euclidean space. In passing from Minkowski 
to Euclidean space we make the substitutions: 

A i A — "l" O O 

d XM = d, XE, i'tpM = V'_B' Pm = -Pe^ 

70M — liEi liM — iliE, IbM = IbE, 

Cm ^ 4'{i^ - MU'"-)il) ^ CE^i)\i^+iMW)il}. (20) 

We write the partition function to which the full QCD partition function 
is rpdiirpH at low momenta as a functional integral over both quark and pion 
fields: B§ 

Z = I Dtt"^ I Dijj'^Dtl^eyi-p ld'^x[ij\x)i^ip{x) 

•d^hd^k2 



(27r)8 



e'^'^'-'^^'^'^ i^Hki)y^M{k^U^^{x)^/M{k2)^p{k2) .(21) 



Eq. (|2l| ) shows quarks interacting with chiral fields U{x), with formfactor func- 
tions equal to the square root of the dynamical quark mass attributed to each 
vertex where U{x) applies. If'^ is a Nf x Nf matrix in flavor and a 4 x 4 
matrix in Dirac indices; it is a unity Nc x Nc matrix in color. 

3.5 Digression: chiral symmetry breaking by instantons 

The partition function (|2l| ) has been actually derived in Ref.(12,13) from con- 
sidering light quarks in the QCD instanton vacuum.Ej Instantons - large topo- 
logical fluctuations of the gluon field - proviiLj; a neat and phenomenologically 
successful rnicroscopic mechanism of SCSB.Ej Because of the famous 't Hooft 
zero modes EEI the random instanton ensemble can be viewed as a collection 
of 'impurities' each binding a quark at exactly zero 'energy'. Owing to the 
quantum-mechanical overlap of quark wave functions quarks can 'hop' from 
one instanton to another, and get delocalized. The density of quarks at zero 
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'energy' becomes finite; this_density is exactly the chiral condensate {'ipip) via 
the Banks-Casher formula.EZI In more physical terms, when 'hopping' from in- 
stanton to anti-instanton and so forth, quarks change their helicity or chirality; 
delocalization means that an infinite number of such flippings are involved, and 
that generates the 'non-slash' term in the quark propagator, i.e. the dynamical 
mass M{k). This quantity is directly related to the Fourier transform of the 't 
Hooft zero mode in the field of an instanton. One can evaluate basic quantities 
related to the SCSB, namely {^ip), M{k) and Ft^ through the characteristics of 
the instanton ensemble, that is their density and average size, these quanJ±t; 
in their turn, being related to Aqcd ( |14| ) from a variational calculation.EJl 

Furthermore, instantons induce (non-local) quark interactions; theipJposon- 
ization leads, at low momenta, to the above partition function (pl|),E^ for a 
review see Ref.(19). However, Eq. ( |2l| ) is probably of a more general nature. 
For example, a similar partition function (but without formfactors) follows 
from the Nambu-Jona-Lasinio model£3 However, contrary to the NJL model, 
instanton-induced interaction preserve all the symmetries of the original QCD: 
they are SULiNf)xSUR{Nf) invariant (in case of TV^ = Nf = 2 they are S'C/ (4) 
invariant) but break explicitly the axial C/4(l) invariance, which is needed to 
solve the C/(l) problem.BEl 

3.6 Local field theory at low momenta 

However realistic is the instanton mechanism of SCSB, actually the low-energy 
partition function (|2l| ) is of a very general nature and in principle one needs 
not refer to instantons to use it. To make it even more general one can add 
to Eq. (pi] ) a Z{k) factor in the kinetic-energy term for quarks and an addi- 
tional formfactor function for pions in the interaction term. However, such 
modifications can be removed by a redefinition of the ip and U fields. 

The formfactor functions y/M{k) for each quark line attached to the chi- 
ral vertex automatically cut off momenta at some characteristic scale. In the 
instanton derivation of Eq. ( |2l| ) this scale is the inverse average size of instan- 
tons, 1/p ~ 600 MeV. In the range of quark momenta k ^ 1/p one can neglect 
this non-locality, and the partition function (^ij) is simplified to a local field 
theory: 

Z= Dtt-^ DiP^DtP exp d^x i^'' {x)[i^+ iMU'"'{x)]i;{x). (22) 



One should remember, however, to cut the quark loop integrals at fc « 1/p « 
600 MeV. Notice that there is no kinetic energy term for pions: it appears 
after one integrates over the quark loop, see below. Summation over color is 
assumed in the exponent of Eq. (P^ ) . 
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Eq. (p^) defines a simple local field theory though it is still a highly non- 
trivial one. Its main properties will be established in the next section. We 
shall see that it reproduces the known properties of the E^L discussed above. 
After that we shall move to building the nucleon from Eq. (|2^) |^. 

4 Properties of the effective chiral lagrangian 

Properties of effective theories of quarks interacting with various meson fields 
have been studied by several authors irLthe 80's, most notably by Volkov and 
Eberta and Dhar, Shankar and Wadia.Q The fact that integrating over quarks 
one gets, in particular, the Wess-Zumino term has been first established by 
Diakonov and EidesP see also below. 

Integrating over the quark fields in Eq. (^2|) one gets an effective chiral 
lagrangian (ExL): 

5'eff[7r] ^ -TVclndet {i^^iMV^) . (26) 

The Dirac operator in Eq. ( pq ) is not hermitean: 

D = i^ + iMU'<'% D'^ = i^ - iMU'^'-'^ , (27) 

therefore the effective action has an imaginary part. The real part can be 
defined as 

Re5'eff[7r] = — -^Indct — , 

D'<D = -d^ + M^ -Mi^W^"), DlDo = -d^ + M^. (28) 



"^We have been askeii-sibout the relation of this low-energy theory with that suggested by 
Manohar and GeorgiE3. One can redefine the quark fields 

^ ^ ^' = exp(J7r'^A'^75/2)V', ip^ -^ ii'^' = i'"' exp{in'^X'^-ys/2), (23) 

and rewrite the lagrangian in ( |22[ ) as 

C = ^1''(i^+ Y+^5+ iM)tP' (24) 

with 

V^='-{i^^e+^''^^.0, A^='-{^d^ii'' -(''d^O, C = exp(i7r'4A-^/2) = l/V^ (25) 

which resembles closely the effective lagrangian of Manohar and Gcorgi (the effective chiral 
lagrangian in a similar form has been independently suggested in Ref.(7)). 
The crucial diff'erence is that Manohar and Georgi have added an explicit kinetic energy 
term _F^Tr (9pC/t9,,f7)/4 on top of Eq. (|24h. This seems to be unnecessary as the kinetic 
energy term arises from quark loops, see tHe next section. 
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In the next two subsections we establish the properties of the real and 
imaginary parts of the E^L separately, following Ref. (24) . 

4-.1 Derivative expansion and interpolation formula 

There is no general expression for the functional (|26) fo r arbitrary pion fields. 

2q) can be estimated nu- 



For certain pion fields the functional determinant 
merically, see section 5. However, one can make a systematic expansion of the 
ExL in increasing powers of the derivatives of the pion field, dU. It is called 
long wave-length or derivative expansion. In fact, one can do it better and 
expand the real part of the E^L in powers of 

'^' iU-1) (29) 



p2 + M2 

where p is the characteristic momentum of the pion field. This quantity be- 
comes small in three limiting cases: (i) small pion fields, n"^{x) ^ 1, with 
arbitrary momenta, (ii) arbitrary pion fields but with small gradients or mo- 
menta, p <C M, (Hi) arbitrary pion fields and large momenta, p 2> M. We 
see thus that expanding the E^L in this parameter one gets accurate results 
in three corners of thp_Hilbert space of pion fields. For that reason we call it 
interpolation formulaEj Our experience is that its numerical accuracy is quite 
good for more or less arbitrary pion fields, even if one uses only the first term 
of the expansion in (|2^) , see below. 

The starting point for both expansions is the following formal manipulation 
with the real part of the ExL (pSl). The first move is to use the well-known for- 
mula. In det [operator] = Spin [operator], where Sp denotes a functional trace. 
One can write: 

ReS'cff[7r] = -^ In dot [l - {-d^ + M^)'^M{^U""')] 

= -^Spln [1 - {-d^ + M^)-^M{^W'-)\ 

= _^ jd^x f-^ e-^^-^Tr In [l - {-d^ + M^yHli^V')] e^'^'" 

= —^ / ^^^ / T^Tr ln[l - (fc2 + m2 -2ik-d- O^y^Mi^U''')}!. (30) 

In going from the second to the third line we have written down explicitly 
what does the functional trace Sp mean: take matrix elements of the operator 
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involved in a complete basis (here: plane waves, exjp{ik ■ x)) , sum over all states 
(here: integrate over d'^k/ {2Tr)*) and take the trace in x. 'Tr ' stands for taking 
not a functional but a usual matrix trace, in our case both in flavor and Dirac 
bispinor indices. In going from the third to the last line we have dragged the 
factor exp(ifc • x) through the operator, thus shifting all differential operators 
d ^ d + ik. We have put a unity at the end of the equation to stress that the 
operator is acting on unity, in particular, it does not differenticitfi it. Above is 
a standard procedure for dealing with functional determinants. E3 

The last line in Eq. (pW can be now expanded in powers of the derivatives 
of the pion field: it arises from expanding ( |30| ) in powers of ^U"^'' and of 
2ik ■ d + d^. The first non-zero term has two derivatives, 

ReSi^jM^^fd^xf^Trf^^^ 



(31) 



off L-j 4 J "- - J (27r)4 ^' \k^ + AP 

df^k M' 



= lJd'xTr{d,U^d,U)-m^J^ 



)4 (fc2+M2)2- 



It is the kinetic energy term for the pion field or, better to say, the Weinberg 
chiral lagrangian: atually it contains all powers of the pion field if one substi- 
tutes U{x) = exp(J7r^(a;)A^/F^), see Eq. (g). The proportionality coefficient 
(the last factor in Eq. (|3l|)) is called F^, experimentally, F^ « MMeV. The 
last factor in Eq. (pTI) is logarithmically divergent; to make it meaningful we 



have to recall that we have actually simplified the theory as given by Eq. (21) 
when writing it in the local form (p2). Actually, the dynamical quark mass 
M is momentum-dependent; it cuts the logarithimically divergent integral at 
k « 1/p. Using the numerical values of ja,« 600 MeV and M « 350 MeV 
according to what follows from instantonsO we find 



, f d^k IvP Nc n I , .■, , . 

^" ^ ^^^ y (^ wtmw " ^^^ ^^ wp - (i«« ^^^^) (^2) 

being not a bad approximation to the experimental value of -FV- One can 
at this point adjust the ultra-violet cutoff 1/p to fit exactly the experimental 
value of F^. Actually, the two-derivative term is the only divergent quantity 
in the E^L: higher derivative terms are all finite. 

A more standard way to present the two-derivative term is by using her- 
mitean Nf x Nf matrices L^ — iUd^U'^ . One can rewrite Eq. (Bl]) as 

Re S^^ [it] = ^ jd^x Tr L^L^, L^ = lUd^UK (33) 
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The next, four-derivative term in the expansion of ReS'ofi is (note that the 
metric is Euchdean) 

Re 5^^^ M = -j^ jd^x [2 Tr {d^L,,f + Tr L,,L„L^L,] . (34) 

These terms describe, in particular, the d-wave tttt scattering lengths, and 
other observables. They can be compared with the appropriate terms of the 
general Gasser-Leutwyler expansion, see Eq. ([lO|). In notations of that equa- 
tion we obtain 

L2 = ^ = 1.58-10"^, vs. (1.35±0.3)- 10"^ (exper.) 

2X2 + -^3 — 0, (experimentally compatible with 0), 

being compatible with the experimental data. The next-to-next-jlDplfiading six 
derivative terms following from Eq. (30) have been computed in£jcJ however, 
a detailed comparison with phenomenology is still lacking here. 

We would like to mention an interesting paper t3 where the derivative 
expansion of the E^L has been obtained from the Lovelace-Shapiro dual reso- 
nance model. For reasons not fully appreciated this model for the tttt scattering 
gives coefficients in front of the 4-derivative terms numerically close to those 
following from Eq. (B3), but there is a discrepancy at the 6-derivative level. 

We now turn to the interpolation formula promised in the beginning of 
this subsection. One can start from the last line in Eq. (pG) and expand it in 
powers of M{^U'^'^). It is clear that the. actual expansion parameter will be 
(p9|). In the first non-zero order we geto 



Re5ir--M = t/^'^/(^T^ 



■4Z 

M(^C/T^) 



{k + idf + M2 



1 



(fc + i9)2 + M2 



Mi^W^'') 



(35) 



It will be convenient now to pass to the Fourier transform of the U{x) field 
understood as a matrix, 



C/(p)= d^e'P-'^iUix)-!]. (36) 



The partial derivatives appearing in Eq. ( p5| ) act on the exponents of the 
Fourier transforms of U, W and become corresponding momenta. As a result 
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we get 



Re5rr^°'M = l|_g_/Tr[C/t(p)t/(p)] 



^ / (27r)4 [(k - f )2 + M2] [(fc + f )2 + M2] ■ ^ ' 



At p -^ the last factor becomes F^ (of. Eq. (|3|)), and Eq. @ is nothing 
but the first term in the derivative expansion, Eq. (pT|). However, Eq. (37) 
also describes correctly the functional S'cff[7r] for rapidly varying pion fields 
(with momenta p ^ M) and for small pion fields of any momenta, when one 
can anyhow expand Eq. ( |3Cl| ) in terms of 7r^(a;) and hence in U{x) — 1. The 
logarithmically divergent loop integral in Eq. (p^ should be regularized, as in 

Similarly, one can get the next term in the 'interpolation' expansion which 



will be quartic in U{p), however our experience tells us that already Eq. (37) 
gives a good approximation to the E^L for most pion fields. 

4-2 The Wess-Zumino term and the baryon number 

We now consider the imaginary partuaf S'cff [tt] . The first non-zexo term in the 
derivative expansion of ImS'eff[7r] isB the Wess-Zumino term. El It cannot be 
written as a d = 4 integral over a local expression made of the unitary U{x) 
matrices, however the variation of the Wess-Zumino term is local. For this 
reason let us ccpsider the variation of Im S'cff [tt] in respect to the pion matrix 
U{x). We haveEd 

Jim S'cff kl = -iVc^ImlndetD^ ^— ^Spf — (51) rSD^ 

■iN 
= ^Sp[iD^D)-^D'<SD-{DD'')-^D^6D^. (38) 

Now one can put in explicit expressions for D,D^ from Eqs.(p^, Eq). The aim 
of this excercise is to get ^U in the denominators so that an expansion in this 
quantity similar to that of the previous subsection could be used. 
Using the Dirac algebra relations (in Euclidean space) 

{Jtj.,lu} = 2(5^^, ll^l^,, 75 = 71727374 = 7^ 

{75, 7p} = 0, 7^ = 75, Tr (757^7^7^75) = 4e„^^5, (39) 

one gets after expanding Eq. (|3q ) in powers of ^W^ the first non-zero term 
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SImSos[n] = ^^ fd'^x e^f3^sTr {d^U^dpUd^U'^dsU U''SU) . (40) 

It can be easily checked that this expression coincides with the variation of the 
Wess-Zumino term written in the form a d = 5 integral Q 

Im Scs[Tr] 

= ^^ jd^x e„0^5, Tr (u^ d^u) (u^dpU) (u^d^u) (u^dsU) (u^d^u) 

+ higher derivative terms. (41) 



The integrand in Eq. (ttlh is a full derivative, however, to write it explicitly 
one would need some parametrization of the unitary matrix U . The expansion 
of Eq. ( pl| ) starts from the fifth power of it^{x) (it is non-zero only if Nf > 3, 
see subsection 3.2). It is important that, similar to ReS'off, the imaginary part 
is also an infinite series in the derivatives. 

The E^L ( p6| ) is invariant under vector flavor-singlet transformations. 
Therefore, there should be a corresponding conserved Noether baryon cur- 
rent, B^. This current is associated with the imaginary part of Scs only; 
since Im ScR is an infinite series in the derivatives so is the associated Noether 
current i?^. For the Wess-Zumino term ( [4l[) the corresponding charge isu 

B ^ -^^ /d^xe^TcTr (u'' diU) (u^djU) (u^dkU) + higher derivative terms. 

(42) 
The explicitly written term is the winding number of the field U{x). Let us 
briefly explain this notion. 

If TT^{x) — > at spatial infinity so that U{x) ^ 1 in all directions, one can 
say that the spatial infinity is just one point. Eq. ( [f2|) gives then the winding 
number for the mapping of the three-dimensional sphere S^ (to which the fiat 
d ~ S space is topologically equivalent when cxi is one point) to the parameter 
space of the SU{Nf) group. In case Nf — 2 the parameter space is also S^ 
so that the mapping is S*^ i— > 5*^. The topologically non-equivalent mappings 
[/(x), i.e. those which can not be continuously deformed one to another, are 
classified by their winding number, an integer analytically given by Eq. (E4). 
In case of Nf > 2 mathematicians prove that mappings are also classified by 
integers given by the same Eq. (p^). 

There exists a prejudice that the baryon number carried by quarks in the 
external pion field coincides with the winding number of that field. Generally 
speaking it is not so because of the higher derivative terms omitted in Eq. (p2|). 
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Only if the pion field is spatially large and slowly varying so that one can neglect 
the higher derivative terms in Eq. (E2|) one can say that the two coincide. 
Otherwise, for arbitrary pion fields, the baryon number is not related to the 
winding number: the former may be zero when the latter is unity, and vice 
versa. 

To see what is going on here, let us calculate directly the baryoiL_num- 
ber carried by quarks in an external time- independent pion field [/(x)o The 
definition of the baryon charge operator in the Minkowski space is 

B^^Jd^^^^oi^. (43) 

Passing to Euclidean space (which we prefer to work with since functional 
integrals are more readily defined in Euclidean) one has to make a substitution 

■0 — > —iip\ 70 — >■ 74, so that 



B = -— d'^i,^-/iij. (44) 

The baryon charge in the path integral formulation of the theory given by 
Eq. ( p^ ) is then 

B = {B) = --^ /d^x (V'^74?/') = -i fd^xTT{x\-/iii^ + iMU''"-)-^\x) 
= ^ Jd^^ Tr (.„ x| ^^^ ^ ^^^J^ ^ ^^^^^^^ |.„ X) 
d3xr^Tr(x| ' " 



— 00 



2tt uj + iH 

Sp 0{-H) = number of levels with E < 0. (45) 



Here 



H = j4lkdk + M-fiV'- (46) 

is the Dirac hamiltonian in the external time-independent pion field C/(x) and 
is a step function. 

Eq. (|4g) is divergent since it sums up the baryon charge of the whole 
negative-energy Dirac continuum. This divergence can be avoided by subtract- 
ing the baryon charge of the free Dirac sea, i.e. with the pion field switched 
out, Hq = -ii-ikdk + M74: 

B^-i fd'^ / ^ Tr (x|— ^ - ^T— |x) ^ Sp [ei-H) Oi-Ho)] . 

(47) 
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In performing the integration over lo we have closed the lo integration contour 
in the upper semiplane. Had we closed it in the lower semiplane we would 
obtain -Sp[6'(iJ) - 6'(iJo)] which is the same resuU since Sp[6'(if ) + B{-II) - 
eiHo) - d{-Ho)] = 0. 

We have thus obtained a most natural result: the baryon charge of quarks 
in the external pion field is the number of negative-energy levels of the hamil- 
tonian ( p^ (the number of the levels of the free hamiltonian subtracted) . 

One can perform the gradient expansion for the baryon number similarly 
to that of the real part of the E^L. To that end let us write 

B = - d'^ —Tr{x\-^—-^-^-^\x) (48) 



_^ 2tt ' 'uj^ + H^ w^ + H^ 



where 



H^ = -dl + Af 2 - Af 74 {dk U^" ) , Hl = -dl + M^ . (49) 

Calculating the matrix element in the plane-wave basis one gets 



dcj r d^k 
2^ J (27r)3 



B = - Id'^l - I 7,^T^14 



7 • (9 + ik) + MU^^ 



7 • (a + ik) 



-(a + ^k)2 + Af2J ^- ^^^^ 

For slowly varying fields U{x) Eq. {pm can be expanded in powers of dlf^ 
and d (applied ultimately to t/^^). Because of the Tr75... the first non-zero 
contribution arises from expanding the denominator in Eq. (pOl) to the third 
power of 7 • (dW^^). Integrals over uj and k are performed explicitly. After 
some simple algebra one gets 

S = --— J d^yie,jkTr{d^U^djUdkU^U) + higher derivative terms (51) 
247r J 

coinciding with Eq. (^) derived, from the Noether current corresponding to 
the Wess-Zumino term (|l|)ExE3 

It should be stressed that the baryon number carried by quarks in the 
background pion field is equal to the topological winding number of the field 
only if it is a slowly varying one. Mathematically, the reason is that for slowly 
varying fields one can neglect the higher derivative terms in Eqs.([42|, |5l| ). Phys- 
ically, the reason is the foUowingO Imagine we start from a pion field U{x) 
whose winding number is one but whose spatial size is tending to zero. Such a 
field has no impact on the spectrum of the Dirac hamiltonian (Eq) : it remains 
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the same as that of the free hamiltonian, namely it has the upper {E > M) 
and lower {E < —M) Dirac continua separated by the mass gap of 2M . The 
baryon number is zero. 

We now (adiabatically) increase the spatial size of the pion field preserving 
its winding number equal to unity. Since the winding number is dimensionless 
this can always be done. At certain critical spatial size the potential well for 
quarks formed by the external pion field is wide enough so that a bound-state 
level emerges from the upper continuum. With the increase of the width of 
the potential well the bound-state level goes down towards the lower Dirac 
continuum. Asymptotically, as one blows up the spatial size of the pion field 
(always remaining in the winding number equal unity sector) the bound-state 
level travels all the way through the mass gap separating the two continua 
and joins the lower Dirac sea - this is a theorem proven in Ref.(24). At this 
point one would discover that there is an extra state close to the lower Dirac 
continuum (as compared to the free, that is no- field case). Therefore, one would 
say that the baryon number is now unity, - in correspondence to Eqs.(^, ^ ). 

In a general case, however, the baryon number of the quark system is the 
number of eigenstates of the Dirac hamiltonian ( ^6| ) one bothers to fill in. The 
role of the winding number of the background pion field is only to guarantee 
that, if the spatial size of the field is large enough, the additional bound-state 
level emerging from the upper continuum is a deep one: asymptotically it goes 
all the way to the lower continuum. 

5 The nucleon 

All constituent quark models start from assuming that the nearly massless light 
quarks of the QCD lagrangian obtain a non-zero dynamical quark mass M « 
350 MeV. This is due to the spontaneous chiral symmetry breaking (SCSB), 
its microscopic driving force being, to our belief, instantons, see subsection 
3.5. Even if one docs not believe in instantons as the microscopic mechanism 
of SCSB one has anyhow to admit that once light quarks get a dynamical or 
constituent mass M{k) such quarks inevitably have to interact with Goldstone 
pions, as explained in subsection 3.4. 

How strong, numerically, is this interaction? Expanding V^ ( |l9| ) to the 
first order of the pion field we get from Eq. ( p^ ) that the linear coupling of 
pions to constituent quarks is, numerically, quite large: 

M , , 

9^m = Tr-4. (52) 

20 



We would like to emphasize that this is a model-independent consequence of 
saying that quarks get a constituent mass. 

At distances between quarks of the order of 0.5 fm typical for interquark 
separations inside nucleons, neither the one-gluon exchange nor the supposed 
linear potential are as large as the chiral forces. Therefore, it is worthwhile 
to investigate whether the chiral forces alone are able to bind the constituent 
quarks inside nucleons. 

5.1 A note on confinement 

The generally accepted confinement scenario is, verbally, very simple. In the 
pure glue world (without light quarks) it is expected that the Coulomb po- 
tential at small separations between static probe quarks is replaced at large 
separations by a linear potential growing ad infinitum,. A constant 14 ton force 
is supposed to be acting between probe quarks even if one is on the Earth and 
the other one is on the Moon. This is a very radical departure from forces 
decaying with distance, which have been known to physics so far. If proven 
correct, it would be one of the most remarkable discoveries in physics of all 
times. 

There have been several proposals of microscopic scenarios which could 
lead to the linear potential, such as the condensation of magnetic monopoles 
or a vacuum filled by random Z{Nc) vortices. We are not in a position to 
discuss them here but would like to note that, despite enormous efforts in the 
last 25 years, none of them have been worked out mathematically so far, at least 
at the semi-quantitative level the instanton vacuum has been constructed. The 
only theoretical model we know of where one gets a linear potential in more 
than two dimensions is the old example by A.PolyakovO of the d — 2 + 1 
Georgi-Glashow model, but there the non-Abelian color group is broken down 
to U{1) by the Higgs condensate. Therefore, our belief that the asymptotic 
linear potential takes place in pure gluodynamics is grounded not on theory 
but mainly on the results of lattice simulations. 

When one looks carefully into how lattice practitioners come to the con- 
clusion about the linear potential one finids several disturbing details. Some 
of them are discussed in our recent papeiuJ. We mention just one here. The 
standard way to extract the quark potential on a lattice is by measuring large 
Wilson loops W{r, t). One side of the loop is called "r" (separation), the other 
one is "t" (time). If one wants to measure the potential V(r) at separation 
r one needs to take t ^ r, otherwise the measurement will be contaminated 
by many excited states. However, such measurements are practically impos- 
sible in the interesting range of r > 1 fm because the average Wilson loop 
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Plateaus for the effective potential 
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Figure 1: Effective potential Vaf{{r,t) = ln[W{r,t)/W(r,t + a)] as function of t at dif- 
ferent values of r from Ref. (35). The data points correspond (from bottom to top) to 
r = 0.65, 0.98, 1.30, 1.62, 1.95 and 2.25 fm. The inverse coupUng used /3 = 2.635 corre- 
sponds to the lattice spacing a = 0.0541 fm. Courtesy G.Bali. 



becomes so small there that measuring it would require statislir;s many orders 
of magnitude exceeding that available with today computers.Ell 

Therefore, this route is abandoned. Instead, Wilson loops are measured in 
the opposite limit, t ^r. The hope is that the slope —d\n.W{r,t)/dt reaches 
its asymptotic (i -^ cxd) value V{r) already at small values of t. One tries to 
assist this early convergence by taking 'smeared' links along the r sides of the 
loop, which are expected to have better overlap with the ground state oL±he 
potential. One of the best (if still not the best) measurement of that kindEj is 
presented in Fig. 1. 

One can see from Fig. 1 that the quality of the plateaus at r < 1 fm 
is quite good though the data tend to slope down as t increases. At larger 
r this trend becomes more pronounced, until the error bars explode so that 
one can hardly get to any conclusions at all. The procedure of extracting the 
potential is detailed in Ref. (35) but basically it follows from the data in Fig. 
1. Thus, the extraction of the linear potential at large r fully relies on the 
faith in that the 'plateaus' of Fig. 1 are not going to slope down as t increases 
from the region t <ti r where they are measured to t ^ r where they should be 
measured. 

As explained inEj this reliance is potentially dangerous. At t ^ r one 
can view the Wilson loop's long side r as 'time' and the short side t as 'sep- 
aration'; then the linear dependence of \nW{r,t) on the long side r, i.e. the 
'linear potential' is built in, whereas it is exactly what is so demanding to 
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prove. This remark becomes even more disturbing in view of the fact that the 
hnear potential is also observed in two cases where the potential has to flatten 
owing to screening: (i) adjoint sources and (ii) fundamental sources but with 
light dynamical fermions. In all cases the same procedure of extracting the 
potentials has been used, namely from measuring Wilson loops at t <ti r, and 
no clear signals of flattening has been observed. 

In our mind, the problem of demonstrating the asymptotic linear potential 
both theoretically and 'experimentally' (on the lattice) remains as challenging 
as it was 25 years ago. 

Whatever the outcome of those very important studies, the world we live in 
is not pure glue: we have light quarks. Light quark-antiquark pairs produced 
from the vacuum eventually screen any rising potential. In the chiral limit 
pions are massless, hence the screening starts at the point where the potential 
exceeds zero energy. In reality the potential cannot exceed 2to^ = 280 MeV 
which is reached, in the standard (Coulomb + linear potential) parametrization 
at a rather small 0.3 fm separation. Nucleon is a stationary state; there is an 
infinite time for any string exceeding that length to decay into pions. Therefore, 
it seems to be senseless to employ large confining potentials in a realistic model 
of nucleons. 

Why, then, do not quarks 'get away' but exist only in bound states? A sim- 
ple explanation can be provided by the SCSB itself. As a result of SCSB origi- 
nally massless or nearly massless quarks get a dynamical mass M « 350 MeV, 
and pions become the lightest excitations in the spectrum. Therefore, light 
quarks need not exist as asymptotic states: instead of producing a quark- 
antiquark pair it is energetically favorable to produce one or several pions. 
Mathematically, it would correspond to the quark propagator with momentum- 
dependent mass, having singularities on the second Riemann sheet under the 
cut starting from the pion threshold. 

What about confinement of heavy quarks (c, 5, t)? Let us imagine, for the 
sake of an argument, that the static potential between probe quark sources in 
a pure glue theory is not rising to infinity but levels off at some V^oo — 2AMoo. 
If AMoo happens to be larger than approximately the light constituent quark 
mass M n heavy quarks would be unstable under decay to D or B mesons. 
Heavy quarks can be thus confined too. The case is to some extent similar 
to electrodynamics with charges Z > 137: such particles are unstable under a 
spontaneous production of e"'"e~ pairs and therefore cannot exist as asymptotic 
states. I 

''This condition is readily met e.g. in the instanton vacuumO rn 

'^The "Z > 137" scenario of confinement has been advocated by V.N.GribovE3 
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Figure 2: Quarks in a nucloon interacting via pion fields. 



5.2 A qualitative picture of the nucleon 

Leaving aside those intriguing problems, we concentrate on the lowest state 
with baryon number one, i.e. the nucleon. As mentioned above the interquark 
separations in the ground-state nucleon are moderate (order of 0.5 fm) and it 
is worthwhile asking whether the simple E^L (p^ ) is capable of explaining the 
basic properties of the ground-state nucleon. Notice that the expected typical 
momenta of quarks inside the nucleon are of the order of M w 350 MeV, that 
is in the domain of applicability of the low-momentum effective theory (p2|), 
see subsection 3.6. 

The chiral interactions of constituent quarks in the 3-quark nucleon, as 
induced by the effective theory (P^), are schematically shown in Fig. 2, where 
quarks are denoted by solid and pions by dashed lines. Notice that, since 
there is no tree- level kinetic energy for pions in Eq. (p2|), the pion propagates 
only through quark loops. Quark loops induce also many-quark interactions 
indicated in Fig. 2 as well. We see that the emerging picture is, unfortunately, 
rather far from a simple one-pion exchange between the constituent quarks: 
the non-linear effects in the pion field are not at all suppressed. 

At this point one may wonder: isn't the resulting theory as complicated 
as the original QCD itself? The answer is no, the effective low-energy theory 
is a simplification as compared to the original quark-gluon theory, because it 
deals with adequate degrees of freedom. Let us imagine that one would like to 
describe 'low-energy' properties of solid states, for example superconductivity. 
Would working with the underlying theory (QED) be helpful? Not too much. 
We know that the microscopic theory leads, under certain conditions, to the 
rearrangement of atoms into a lattice, so that translational symmetry is spon- 
taneously broken. As a result the Goldstone bosons appear (here: phonons), 
and electrons get a dynamical mass different from the input one. The most 
important forces are due to phonon exchange between electrons: in fact they 
are driving superconductivity in the BCS theory. Playing with this analogy, 
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nucleon is like a polaron (a bound state of electrons in the phonon field) , rather 
than a positronium state in the vacuum. After chiral symmetry is broken we 
deal with a 'metal' phase rather than with the vacuum one, and one has to use 
adequate degrees of freedom to face this new situation. 

Instantons play the role of the bridge between the microscopic theory 
(QCD) and the low-energy theory where one neglects all degrees of freedom ex- 
cept the Goldstone bosons and fermions with the dynamically-generated mass. 
Instantons do the most difficult part of the job: they explain why atoms in 
metals are arranged into a lattice, what is the effective mass of the electron 
and what is the strength of the electron-phonon interactions. However, one 
can take an agnostic stand and say: I don't care how 350 MeV is obtained from 
the microscopic Aqcd and why do atoms form a lattice in the metals: I just 
know it happens. To such a person we would advise to take the low-energy 
theory (^) at face value and proceed to the nucleon. 

A considerable technical simplification is achieved in the limit of large Nc- 
For Nc colors the number of constituent quarks in a baryon is Nc and all 
quark loop contributions are also proportional to Nc, see section 4. Therefore, 
at large Nc one can speak about a classical self-consistent pion field inside 
the nucleon: quantum fluctuations about the classical self-consistent field are 
suppressed as l/Nc- The problem of summing up all diagrams of the type 
shown in Fig. 2 is thus reduced to finding a classical pion field pulling Nc 
massive quarks together to form a bound state. 

5.3 Nucleon mass: a functional of the pion field 

Let us imagine that there is a classical time-independent pion field which is 
strong and spatially wide enough to make a bound-state level of the Dirac 
hamiltonian (|4q ) for massive quarks, call its energy -Eicvoi- We fill in this level 
by Nc quarks in the antisymmetric state in color, thus obtaining a baryon 
number one state, as compared to the vacuum. The interactions with the 
background chiral field are color-blind, so one can put Nc quarks on the same 
level; the fact that one has to put them in an antisymmetric state in color, i.e. 
in a color-singlet state, follows from Fermi statistics. 

One has to pay for the creation of this trial pion field, however. Call this 
energy -Efioid- Since there are no direct terms depending on the pion field in the 
low-momentum theory (E2h the only origin of i^ficid is the fermion determinant 
( P6| ) which should be calculated for time- independent field U{x). It can be 
worked out with a slight modification of the equations of section 4. We have:E3 

Scs[tt] = -A^clndetf — 
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= -Nc Sp [\n[idt + iH) - \n{idt + iH^)] 

/I 
-^ Sp [ln(tj + iH) - ln(w + iHo)] , (53) 

27r 

where H is the Dirac hamiltonian ( jiq ) in the stationary pion field, Hq is the 
free hamiltonian and T is the (infinite) time of observation. Using an important 
relation ^ 

SpiH -Ho) = (54) 

(telling us that the sum of all energies, with their signs, of the Dirac hamil- 
tonian (Bq) is the same as for the free case) one can integrate in Eq. (m) by 
parts and get 



S'off[7r] = TEficid = TNc --Sp 

J ^TT 



uj + iH uj + iHo 
= TN, J2 [En-E^^A. (55) 



Going from the first to the second line we have closed the to integration contour 
in the upper semiplane; owing to the trace relation (p3) closing it in the lower 
semiplane would produce the same result. 

We see that the price iJficid one pays for a creation of the time- independent 
pion field coincides with the aggregate energy of the lower Dirac continuum 
in that field. The energy of the additional level emerging from the upper 
continuum, which one has to fill in to get the baryon number one state ,|_Ei_,oij 
should be added to get the total nucleon mass. This simple scheme EI1l3 is 
depicted in Fig. 3. Naturally, the mass of the nucleon should be counted 
from the vacuum state corresponding to the filled levels of the free lower Dirac 
continuum. Therefore, the (divergent) aggregate energy of the free continuum 
should be subtracted, as in Eq. (|55|). 

We have thus for the nucleon mass: 

Mn= min (iVc£;icvei[7r]+£;ficidM). (56) 

Both quantities, -Biovci and -Eficid, are functionals of the trial pion field 
TT (x). The classical self-consistent pion field is obtained from minimizing the 

''It follows from taking the matrix trace of the hamiltonian (llQ). It should be kept in mind, 
though, that such a naive derivation can be potentially dangerous because of anomalies in 
infinite sums over levels. However, it can be checked that in this particular case there are 
no anomalies and the naive derivation is correct. 
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Figure 3: Spectrum of the Dirac hamiltonian in trial pion field. The solid lines show occupied 
levels. 



nucleon mass (|6|) in 7r^(x). It is called the soliton of the non-linear functional 
(p6|), hence the Chiral Quark-Soliton Model. An accurate derivation of 



Eq. (56) from the path- integral representation for the correlation function of 
loffe currents is presented in Ref.(24). It solves the problem of summing up all 
diagrams of the type shown in Fig. 2 in the large- iVc limit. 

The idea that a sigma model with pious coupled directly to constituent 
quarks, can be used to bi. 



jd the nucleon soliton has been first suggested by 

Kahana, Ripka and SoniEa and independently by Birse and Banerjee.E3 We 
would call them the authors of the Chiral Quark-Soliton Model. Technically, 
however, in these references an additional ad hoc kinetic eneEgAj-term for pion 
fields has been used, leading to a vacuum instability paradox.cJ'cll The present 
formulation of the model has been given in Ref.(37), together with the discus- 
sion of its domain of applicability and its physical contents. A detailed theory 
based on the path integral approach paving the way for calculating nucleon 
observables has been presented in Ref.(24). 

5.4 Nucleon profile 

To find the classical pion field minimizing the nucleon mass (^6|) one has first 
of all decide on the symmetry of the pion field. Had the field been a singlet 
one would take a spherically-symmetric ansatz. However, the pion field has 
fiavor indices A = 1, ..., iV? — 1. At Nf ~ 2 the three components of the pion 
field can be married with the three space axes. This is called the hedgehog 
ansatz; it is the minimal generalization of spherical symmetry to incorporate 
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the 7r± = (ttI ± iiT'^)/^2 and tt^ ^ vr^ fields: 

TT^{^)/F, = n^P{r), n^ = — , r = |x|, (57) 

r 

where P{r) is called the soliton profile. 

The choice of the ansatz is not at all innocent: baryons corresponding to 
different choices would have qualitatively different properties, see the next sub- 
section. The maximally-symmetric ansatz (p7^ will have definite consequences 
in applications. 

In the Nf — 3 case there are 8 components of the 'pion' field, and there 
are several possibilities to marry—them to the space axes. The commonly used 
is the 'left upper corner' ansatz:E2l 



C/(x) ^ exp (z7r-(x)A-) = | °"P t^^" ' ^^^^'^^ ] , n = ^. (58) 

1 / ^ 

As we shall see in the next subsection, the quantization of rotations for this 
ansatz leads to the correct jSpectrum of the lowest baryons. 

Another SU{3) ansatzcJ discussed in the hterature is (/, g,h — 1,2,3 are 
the flavor indices): 



Ufg = e'^=/3 ^^og p^ g^^ j^ ^g-^P2 _ j^os p^) „^.„^ ^ sjj^ p^ efghUh] , (59) 

where Pi^2{t) are spherically-symmetpc profile functions. This ansatz is used 
to describe strangeness —2 dibaryons.c3 We shall not consider it here but con- 
centrate on the usual baryons for which the hedgehog ansatz (B^) or (BSh is 
appropriate. 

Let us first discuss restrictions on the best profile function P(r) which 
should mimimize the nucleon mass (|5^). What is the asymptotics of P{r) at 
large r? To answer this question one has to know the behavior of iJfieid [tt] for 
slowly varying pion fields. Using Eq. ( p5| ) as a starting point one can work out 
the derivative expansion of the furictipnal Efidd [ti"] similar to that for the full 
ExL, see subsection 3.1. We haveEZllHJ 

p2 f JV f 

£^ficidM = -f d^^TTL,L, - ^^ /d^x [2Tr{d,Lif + TYL,L^L,Lj] 

+ higher derivative terms, Li = iUdiU^ . (60) 
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Substituting here the hedgehog ansatz one gets a functional of the profile 
function P{r); varying it one finds the Euler-Lagrange equation valid for slowly 
varying profiles, in particular, for the tail of P{r) at large r. It follows from 
this equation that P{r) — Ajr^ at large r. The second contribution to the 
nucleon mass, -Eiovci [""] , does not alter this derivation since the bound-state 
wave function has an exponential but not power behavior at large r. Actually 
we get the pion tail inside the nucleon, and the constant A is related to the 
nucleoiL^xial constant. This relation is identical to the one found in the Skyrme 
modehEil 

gA = ^AFl. (61) 

The exponentially decreasing wave function of the bound-state level does not 
change this equation, as well as the Goldberger-Treiman relation for the pion- 
nucleon coupling constant, 

a^NN = —^ = ^ -^- (62) 



Furthermore, it follows from the next four-derivative term in Eq. (pOl) that the 
1/r^ correction to Pir) at large r is absent! OE^ It means probably that the 
pion tail inside the nucleon is unperturbed to rather short distances. 

The second important question is what should we choose for P(0)? As 
explained in subsection 3.2, a quantity which guarantees a deeply-bound state 
in the background pion field is the winding number of the field, Eq. ([42|). 
Substituting the hedgehog ansatz into Eq. (B^) one gets 



2 r...„2..„.rfP(r)^_l 
dr n 



2 f°" 

^wind = / dr sin^ P(r) 

^ Jo 



Pi^r) - ^^"^^^"^ 



(63) 



Since P{oo) = the way to make this quantity unity is to choose P(0) = tt = 
3.141592... 

An example_n£^a one-parameter variational function satisfying the above 
requirements isLllEj 

P{r) = 2 arc tan (—) , A = 2r^. (64) 



The Dirac hamiltonian (Uq) in the hedgehog pion field (|57| ) commutes 
neither with the isospin operator T nor with the total angular momentum 
J = L + S but only with their sum K — T + J called the 'grand spin'. 
The eigenvalue Dirac equations for given value of K^jK^ have been derived 
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in Rcf.(24). Generally speaking, there appears a bound-state level with the 
K^ = 0+ quantum numbers whose energy can be found from solving the 
Dirac equations for two spherically-symmetric functions j, h, 

-M sin Ph+ {Eicvc\ + M cos P) j, 

M sin P j + (-i'lcvci + M cos P) h, (65) 

with the boundary conditions h{0) = 1, j(0) = Cr, h{oo) = j{oo) = 0. 

These equations determine one of the two contributions to the nucleon 
mass, -Eiovci- The second contribution, namely that of the aggregate energy of 
the lower Dirac continuum in the trial pion field, which we have called Encid, 
can be found in several different ways. One way is to find the phase shifts in 
the 1nwpr r-nntimnim arising from solving the Dirac equation for definite grand 
spin i^rEflc^ Another method i^o diagonalizc the Dirac hamiltonian (^) in the 
so-called Kahana-Ripka basiscil written for a finite- volume spherical box. Both 
methods_axe, numerically, rather involved. There exists a third (approximate) 
method EZrcJ allowing one to make an estimate of -Bfidd in a few minutes on a 
PC. It is based on the interpolation formula for the E^L, see Eq. (|37|). 

Let us discuss the qualitative behavior of -Biovci and i^ficM with the soli- 
ton scale parameter rg assuming for definiteness that the profile is given by 
Eq. (H). 

The trial pion field plays the role of the (relativistic) potential well for 
massive quarks. The 'depth' of this potential well is fixed by the condition 
P(0) — TT and is always finite: this is related to the fact that the pion field has 
the meaning of angles. The spatial size of the trial pion field tq plays the role 
of the 'width' of the potential well. It is well known that in three dimensions 
the condition for the appearance of a bound state is MVr'^ > const where V is 
the depth of the well and 'const' is a numerical constant of the order of unity 
depending on its concrete shape. In our case V « M, so the condition that 
the bound state appears is tqM = 0{1). Therefore, at small sizes tq there is 
no bound state for the Dirac hamiltonian (E6f), so that i^icvci coincides with 
the border of the upper Dirac continuum, -Bicvei = +M. At certain critical 
value of To a weakly bound state emerges from the upper continuum. [For the 
concrete ansatz (bih the threshold value is roAI ~ O.5.] As one increases tq 
the bound state goes deeper and E'lcvci monotonously decreases. At very large 
spatial sizes, tq — *■ c», iiwoi approaches the lower continuum, its difference 
from — M falling as l/r^t^ The behavior of i^icvci as function of ro is plotted 
in Fig. 4. 

The monotonous decrease of E'lovci with the increase of ro is a prerogative 
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Figure 4; Nucleon mass and its constituents as function of the soliton size. The short-dash 
line shows S-Bieveh the long-dash line is -Efleldi the solid line is their sum, A^jv- 



of the trial pion field with winding number 1. Had it been zero, -E'lovci would 
first go down and then start to go up, asymptotically joining back the upper 
continuum. In the case of iVwind = — 1 the bound state would travel in the 
opposite direction: from the lower towards the upper continuum. At A'wind = n 
as much as n levels would emerge, one by one, from the upper continuum and 
travel all the way through the mass gap towards the lower one. For the trial 
pion field of the hedgehog form all these things happen exclusively for states 
with grand spin iiT = 0.cJ 

Turning now to i^fioid we first notice that for large spatial sizes of the trial 
pion field one can use the first term in the derivative expansion for -Eficid, see 
Eq. (|6C|). On dimension grounds one immediately concludes that -Eficid ~ -f^^o 
for large rg, i.e. is infinitely linearly rising in rg. At small tq a slightly more 
complicated analysis Ej shows that E'fidd ~ ^o- On the whole, -Eficid is a 
monotonously rising function of tq shown in Fig. 4. 

The nucleon mass, M.m ~ 3£'ievci + -E-fieid (for Nc = 3) is also plotted 
in Fig. 4 taken from Refs. (37,24). One observes a non-trivial minimum for 
Mn corresponding to rp ~ 0.98/Af ~ 0.57 fm. This is, phenomenologically, 
a very reasonable value, since from Eqs.(0, |62[ ) one immediately gets gA — 
1.15 versus 1.25 (exp.) and g-nNN — 13.6 versus 13.5 (exp.). The nucleon mass 
appears to be Mn ^ HOO MeV with 3^iovoi ^ 370 MeV, ^fioid ^ 730 MeV. 
Note that the 'valence' quarks (sitting on the bound-state level) come out to 
be very strongly bound: their wave function falls off as exp(— r/0.6 fm), and 
about 2/3 of the quark mass M ~ 350 MeV is eaten up by interactions with 
the classical pion field. Relativistic effects are thus essential. 
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Though the nuclcon bound state appears to be somewhat higher than the 
free-quark threshold, 3M ~ lOSOMeV, there are several known corrections to it 
which are negative. The largest correction to the nucleon mass is due to taking 
into account explicitly the one-gluon exchange between both 'valence' and 'sea' 
quarks; this correction is 0{Nc) aa^is the nucleon mass itself. Numerically, it 
turns out to be about —200 MeVcZl and seems to move the nucleon mass just 
into the right place. Q 

We see thus that the 'valence' quarks in the nucleon get bound by a self- 
consistent pion field whose energy is given by the aggregate energy of the 
negative Dirac continuum distorted by the presence of the external field. This 
picture of the nucleon interpolates between the old non-relativistic quark mod- 
els (which would correspond to a shallow bound-state level and an undistorted 
negative continuum) and the Skyrme model (which would correspond to a spa- 
tially very large pion soliton so that the bound-state level would get close to 
the lower continuum and the field energy iJficid would be given just by a couple 
of terms in its derivative expansion) . The reality is somewhere in between: the 
bound-state level is a deep one but not as deep as to say that all the physics 
is in the lower Dirac continuum. 

Idelogically, this picture of the nucleon at large Nc is somewhat similar 
to the Thomas-Fermi picture of the atom at large Z. In that case quantum 
fluctuations of the self-consistent electrostatic field binding the electrons are 
suppressed by large Z, however corrections go as powers of Z~^/^, as contrasted 
to the N~^ corrections of the Chiral Quark-Soliton Model. Therefore, the 
latter model is in a slightly better position in respect to quantum corrections 
than the former. 

Apart from using the large Nc approximation (which is in fact just a tech- 
nical device needed to justify the use of the classical pion field) the Chiral 
Quark-Soliton Model makes use of the small algebraic parameter {Mpf « 1/3 
where 1/p is the UV cutoff of the E^L (|2^). [In the instanton derivation p is 
the average size of instantons in the vacuum.] This p is, roughly, the size of 
the constituent quark, while the size of the nucleon is, parametrically, 1/Af . 
The fact that the constituent quark picture works so well in the whole hadron 
physics finds its explanation in this small numerical parameter. [In the instan- 
ton picture it is due to the relative diluteness of the instanton vacuum, which 
in its turn is related-lo the 'accidentally' large number (11/3) in the asymp- 
totic freedom law.EJIla] The small parameter (Mp)^ <C 1 makes it possible to 
use only quarks with dynamically generated mass and chiral fields as the only 



'^There exist also numerous quantum corrections to the nuclcon mass of different origin, 
which are of the order of O(Af^). Unfortunately, it is difficult today to treat them in a 
systematic fashion. 
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essential degrees of freedom in the range of momenta k ^ M ^ 1/p, and that 
is exactly the range of interest in the nucleon binding problem. 

The ahmffi numerics have been obtained from the interpolation formula 
for E^fioid -EZtcJ Exact calculations of iJfioid performed in c3 as well as taking 
more involved profiles with three variational parameters did not lead to any 
significant changes in th e numerics. 

Following refs. EZio there had been many calculations of the nucleon mass 
and of the 'best' profile_iising various regularization schemes and parameters 
of the chiral model, seec3 for a review. The effective low- momenta theory (22) 



comes along with an intrinsic ultraviolet cutoff, in the form of a momentum 
dependence of the constituent quark mass, M{k). In the instanton approach, 
it can be shown on general grounds that this is a rapidly decreasing function 
at momenta of the order of the inverse average instanton size, 1/p ~ 600 MeV. 
However, the present 'state of the art' does not allow one to determine this 
function accurately at all values of momenta - to do so, one would need a very 
detailed understanding of the QCD vacuum. This places certain restrictions 
on the kinds of quantities which can sensibly be computed using the effective 
theory. Those are either finite ones, which do not require an UV cutoff at all, 
or quantities at most logarithmically divergent. Both type of quantities are 
dominated by momenta much smaller than the UV cutoff, k <^ 1/p, so one 
can compute them mimicking the fall-off of M{k) by an external UV cutoff 
A ~ 1/p, using some regularization scheme. V\ Fortunately, almost all nucleon 
observables belong to these two classes. The uncertainty related to the details 
of the ultra-violet regularization leads to a 15-20% numerical uncertainty of 
the results, and that is the expected accuracy of the model today. 

6 Quantum numbers of baryons 

The picture of the nucleon outlined in the previous section is "classical" : the 
quantum fluctuations of the self-consistent pion field binding Nc quarks are 
totally ignored. Among all possible quantum corrections to the nucleon mass 
a special role belongs to the zero modes. Fluctuations of the pion field in the 
direction of the zero modes cannot be considered small, and one has to treat 
them exactly. Zero modes are always related to continuous symmetries of the 
problem at hand. In our case there are 3 zero translational modes and a certain 
number of zero rotational modes. The latter determine the quantum numbers 
of baryons; it is here that the hedgehog (or whatever) symmetry of the ansatz 
taken for the self-consistent pion field becomes crucial. 



•^Recently the sclf-cxmspticnt field binding the nucleon has been found using directly the 
non-local coupling (Elh-H 
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A general statement is that if the chiral field C/ci(x) minimizes the nucleon 
mass functional (pq), a field corresponding to rotated spatial axes, Xi — > OijXj, 
or to a unitary-rotated matrix in flavor space, Ud — > RUdR'' , has obviously 
the same classical mass. This is because the functional (^6|) to be minimized 
is isotropic both in flavor and ordinary spaces. 

Specifically for the hedgehog ansatz [see Eq. ( |57| ) for the flavor SU{2) and 
Eq. (^) for the SU{3)] any spatial rotation is equivalent to a flavor rotation. 
We shall first show it for a more general case of SU{3). 

6.1 Case of three flavors 

The space-rotating 3x3 matrix Oij can be written as 

0,, = ^Ti{SnS^T,) (66) 

where S is an SU{2) 2x2 matrix and r.; are the three Pauli matrices. One 
can immediatelly check that Oij are real orthogonal 3-parameter matrices with 
OijOkj = Sik and OijOik = Sjk, as it should be. 

When one rotates the space putting n'^ = rijOji the 2x2 matrix standing 
in the left upper corner of the ansatz (p8|) can be written as 



exp [i{n' ■ x)P(r)] = cosP(r) + i(n' • t) sinP(r) 

= S [cos P{r) + i{n ■ t) sin P{r)] SK (67) 

Therefore, if one consideres the hedgehog ansatz (^^ rotated both in flavor and 
usual spaces, the latter can be completely absorbed into the former one: 

R [/ci(Ox) i?t = ^ L/ci(x) 7?t (68) 

with 

/ S ° 

R = R\ I . (69) 

\ f 

For that reason it is sufficient to consider rotations only in the flavor space. 
Hence there are 3 zero rotational modes in the SU{2) and 8-1=7 in the SU{3) 
flavor case. The rotation of the form R — exp(iaA*) commutes with the left- 
upper-corner ansatz and therefore does not correspond to any zero mode. This 
will have important consequences in getting the correct spectrum of hyperons. 

The general strategy is to consider a slowly rotating ansatz 
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t/(x,t) = i?(i)(7ci(x)i?t(i) (70) 

and to expand the energy of the bound-state level and of the negative Dirac 
continuum in 'right' {^Ia) and 'left' (17^) angular velocities 



fJA = -«Tr(i?^i?A^), fiA = -iTr(i?i?tA^), n^ ^ n^ ^ 2Ti R'' R. (71) 

Taking into account onlv_th£j lowest powers in the time derivatives of the 
rotation matrix R{t) one getscJ^Ej the following form of the rotation lagrangian: 

i"' = ^lAB^AnB - ^f^8. (72) 

Here Iab is the SU{3) tensor of the moments of inertia, 

/.. = ^/^Trf^A-^A-) (73) 

where the uj integration contour should be drawn above the bound-state energy 
-Eiovci to incorporate the 'valence' quarks. 

The appearance of a linear term in fig is an important consequence of 
the presence of an extra bound-state level emerging from the upper Dirac 
continuum, which fixes the baryon charge to bejinity. In the Skyrme model this 
linear term arises from the Wess-Zumino term.LiI For simplicity we have written 
unregularized moromts of inertia, though Eq. ( [73[ ) should be regularized in 
some way, see e.g. Ej 

Owing to the left-upper-corner ansatz for the static soliton the tensor Iab 
is diagonal and depends on two moments of inertia, /i 2: 



h Sab, 


A,B = 1,2,3, 


h Sab, 


^,5 = 4,5,6,7, 


0, 


A,B = 8. 



Iab^< I26AB, AS = 4,5,6,7, (74) 

[ 0, A,B = 8. 

Therefore, the rotational lagrangian ( [72| ) can be rewritten as 

2 ti 2 ^4 2V3 ' 



To quantize this rotational Lagrangian one 
zation procedure, same as in the Skyrme model 



|ifj3liiiii 



canonical quanti- 
Introducing eight 
angular momenta canonically conjugate to 'right' angular velocities flA, 
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a rrot 

and writing the hainiltonian as 

if™' = VLaJa ~ i™* (77) 

one gets 

^ A=l ^ A=4 

with the additional quantization prescription following from Eq. (uq), 

In the Skyrme model this quantization rule follows from the Wess-Zumino 
term. In our approach it arises from filling in the bound-state level, i.e. from 
the 'valence' quarks. It is known to lead to the selection rule: not all possible 
spin and SU{3) multiplets are allowed as rotational excitations of the SU{2) 
hedgehog. Eq. (JT^) means that only those SU{3) multiplets are allowed which 
contain particles with hypcrchargc y = 1; if the number of particles with 
y = 1 is denoted as 2J + 1, the spin of the allowed SU{3) multiplet is equal 
to J. 

Therefore, the lowest allowed SU{3) multiplets are: 

• octet with spin 1/2 (since there are two baryons in the octet with y = 1, 

the N) 

• decuplet with spin 3/2 (since there are four baryons in the decuplet with 
y = 1, the A) 

• anti-decuplet with spin 1/2 (since there are two baryons in the anti- 
decuplet with y = 1, the iV*) 

The next are 27-plets with spin 1/2 and 3/2 but we do not consider them here. 
We see that the lowest two rotational excitations are exactly the lowest 
baryon multiplets existing in reality. The third predicted multiplet, the anti- 
decuplet, contains exotic baryons which cannot be made of three quarks, most 
notably an exotic Z+ baryon having spin 1/2, isospin and strangeness -1-1. A 
detailed study of the anti-decuplet performed recently ino predicts that such 
a baryon can have a mass as low as 1530 MeV and be very narrow. Several 
experimental searches of this exotic baryon are now under way. 
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It is easy to derive the splittings between the centers of the multiplets 
listed above. For the representation (p, q) of the SU{3) group one has 

8 1 

E'^a = :t[p' + '?'+M + 3(p + 9)], (80) 

A=l "^ 

therefore the eigenvalues of the rotational hamiltonian d7q)) are 

i^(-) = g^[p^ + .^+M + 3b + ,)]+(^-^)j(J + l)-A. (81) 
We have the following three lowest rotational excitations: 



(p, q) = 


= (1,1), 


J =1/2: 


octet, spin 1/2, 


(82) 


(p, q) = 


= (3,0), 


J = 3/2 : 


decuplet, spin 3/2, 


(83) 


(p, q) = 


= (0,3), 


J=l/2: 


anti-decuplet, spin 1/2. 


(84) 



The splittings between the centers of these multiplets are determined by the 
moments of inertia Ii^2- 

3 

Aio_8 = ^(3*0) - ^(1*1) = ^, (85) 

3 

^To-8 = -^(o'a) ^ ^(° 1) = ^- (86) 

The appropriate rotational wave functions describing member s of these 
multiplets are given by Wigner finite-rotation functions D®'^°'^°(i?).E3£J 

When dealing with the flavor SU{3) case neglecting the strange quark 
mass rUs is an over-simplification. In fact, it is easy to incorporate rris 7^ in 
the first order. As a result one gets very reasonable sphttings iriaidnJLiie SU{3) 
multiplets, as well as mass corrections to different observables.E2lc3C3 

In general, the idea that all light baryons are rotational excitations of one 
object, the 'classical' nucleon, leads to numerous relations between properties 
of the members of octet and decuplet, which follow purely from symmetry 
considerations and which are all satisfied up to a few percent in nature. The 
SU{3) symmetry by itself says nothing about the relation between jiifferent 
multiplets. Probably the most spectacular is the Guadagnini formulaEll which 
relates splittings inside the decuplet with those in the octet, 

8(mH* -I- mpf) + 3toe = IIttia + Stoe*, (87) 

which is satisfied with better than one-percent accuracy! 
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6.2 Case of two flavors 

If one is interested in baryons predominantly 'made of u, d quarks, the flavor 
group is SU(2) and the quantization of rotations is more simple. 
In this case the rotational lagrangian is just 



3 3 

-'1 V^^2 -'1 



2 



i=l A=l 

where the 'right' (fli) and 'left' (Qa) angular velocities are defined by Eq. ([7l|). 
This is the lagrangian for the spherical top: the two sets of angular velocities 
have the meaning of those in the 'lab frame' and 'body fixed frame'. The 
quantization of the spherical top is well known from quantum mechanics. One 
has to introduce two sets of angular momenta, Si (canonically conjugate to 
fii) and Ta (conjugate to (Ia)- Both sets of operators act on the coordinates 
of the spherical top, say, the Euler angles. It will be more convenient for us to 
say that the coordinates of the spherical toe are just the entries of the unitary 
matrix R defining its finite-angle rotationtj 

The angular momenta operators Si, Ta act on R as generators of right 
(left) multiplication, 

g»(aS)jr^g-i(aS) ^j^gi(aa)^ (89) 

e«(oT)^g-»(aT) ^ g-^(ar) j^ (gg) 

and satisfy the commutation relations 

[Ta,Tb] = ieABcTc, [Si, Sj] — ie^jkSk, 

[Ta,S,]^0, {TAf = {S,f. (91) 

A realization of these operators is 



Si — R, 



(1l 
\ 2 



d_ 



kq ORpq 

d 



The rotational hamiltonian is 



7J'°* = nA - L'°* = ^aTa ~ U°' = 1^ = ^. (93) 



2/i 21 
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Comparing the definition of the generators ( |89| , |90D with the ansatz (70) we 
see that Ta is the flavor (here: isospin) operator and Si is the spin operator, 
since the former acts to the left from R and the latter acts to the right. 

The normalized eigenfunctions of the mutually commuting operaiors 6*3, T^ 
and S^ = T^ with eigenvalues S3, T3 and S{S + 1) = T{T + 1) areH 



4^^3^)(i?) = V2S + Ti-lf+^-wi%2(^) (94) 

where D{R) are Wigner finite-rotation matrices. For example, in the S = T = 



1 /? 

1/2 representation Dpq (i?) = Rpq, i.e. coincides with the unitary matrix R 



itself. 

The rotational energy is thus 

2/1 2/1 

and is (2S'+ 1)^ = (2r-|- l)^-fold degenerate. The wave functions ( |94| ) describe 
at S' = T = 1/2 four nucleon states (proton, neutron, spin up, spin down) and 
at S* = T = 3/2 the sixteen A-resonance states, the splitting between them 
being 

mA-rriN^^^OiN-') (96) 

All 

(coinciding in fact with the splitting between the centers of decuplet and octet 
in the more general SU{S) case, see Eq. (^). 

It is remarkable that the nucleon and its lowest excitation, the A, fits into 
this spin-equal-isospin scheme, following from the quantization of the hedge- 
hog rotation. Moreover, since N and A are, in this approach, just different 
rotational states of the same object, the 'classical nucleon', there are certain 
relations between their properties. These relations are identical to those found 
first in the Skyrmc model Eil since they follow from symmetry considerations 
only and do not depend on concrete dynamics which is of course different in 
the naive Skyrme model. For example, one getsJor the dynamics-independent 
ratio of magnetic moments and pion couplings EiJ 

^J■AN 1 



Mp ~ Mn V ^ 



0.71 vs. 0.70 ±0.01 (exp.), 



^^^^^ = ^ = 1.5 vs. 1.5 ±0.12 (exp.). (97) 

gyrNN 2 

We should mention that there might be interesting implications of the 
'baryons as rotating solitons' idea to nuclear physics. The low-energy interac- 
tions between nucleons can be viewed as interactions between spherical tops 
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depending on their relative orientation RiRI in the spin-isospin spacesLjEJ It 
leads to an elegant description of NN and A^ A interactions in a unified fashion, 
and it would be very interesting to check its experimental consequences (as far 
as we know this has not been done yet). A nuclear medium is then a medium 
of interacting quantum spherical tops with extremely anisotropic interactions 
depending on relative orientations of the tops both in the spin-isospin and in 
ordinary spaces. i_. 

This unconventional point of view is supported by the observation l3 that 
one can get the correct value of the so-called symmetry energy of the nucleus, 
(25MeV) • {N~Zy/A, the numerical coefficient appearing as l/8/i ~ 25MeV 
where /i is the SU{2) moment of inertia; from the A — N splitting ( pq ) one 
finds /i ~ (200 MeV)~^. We do not know whether the language of spher- 
ical tops is fruitful to describe ordinary nuclear matter (probably it is but 
nobody tried), however it is certainly useful to address new questions, for ex- 
ample whether nuclear mattfir at high densities can be in a strongly correlated 
antiferromagnet-type phase. E3 

Finally, let us ask what the next rotational excitations could be? If one 
restricts oneself to only two flavors, the next state should be a (5/2, 5/2) 
resonance; in the three-flavor case the next rotational excitation is the anti- 
decuplet with spin 1/2, see subsection 6.1. Why do not we have any clear 
signal of the exotic (5/2, 5/2) resonance? The reason is that the angular 
momentum J = 5/2 is numerically comparable to Nc ~ 3. Rotations with J « 
Nc cannot be considered as slow: the centrifugal forcp-jdeform considerably 
the spherically-symmetrical profile of the soliton field;c2[E3 simultaneously at 
J fa Nc the radiation of pions by the rotating body makes the total width of 
the state comparable to its mass.c3t3 In order to survive strong pion radiation 
the rotating chiral solitons with J > Nc have to stretch intp-cigar-like objects; 
such states lie on linear Regge trajectories with the slope c3 a' ~ I/Stt'^F^ « 
1.45 GeV~^. One cannot exclude an intriguing possibility that all large- 
hadron states lying on Regge trajectories are actually rotating chiral solitons. 

The situation might be somewhat different in the three-Havor case. First, 
the rotation is, roughly speaking, distributed among more axes in flavor space, 
hence individual angular velocities are not necessarily as large as when we 
consider the two-flavor case with J — 5/2. Actually, the SU{2) baryons with 
J = 5/2 belong to a very high multiplet from the SU{3) point of view. Second, 
the radiation by the soliton includes now K and t] mesons which are substan- 
tially heavier than pions, and hence such radiation is suppressed. Actually, the 
anti-decuplet seems to have moderate widths e3 and it is worthwhile to search 
for the predicted exotic states. 
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7 Some applications 

There exists by now a rather vast Uterature studying baryon observables in 
the Chiral Quark-Soliton Model. Baryon formfactors (electric, magnetic and 
axial), mass splittings, the nucleon sigma term, magnetic moments, weak de- 
cay constants, tensor charges and many other characteristics of nucleons and 
hyperons have been calculated in the model. We address the reader to an 
extensive reviewed on these matters. 

Here we would like to point out several developments of the Chiral Quark- 
Soliton Model interesting from the theoretical point of view. The list below is, 
of course, very subjective. 

The study of the spin content o£_the nucleon in the model has been pi- 
oneered by Wakamatsu and Yoshiki.E3 They showed that the fraction of the 
nucleon spin carried by the spin of quarks is about 50% (and could be made 
less): the rest is carried by the interquark orbital moment, the Dirac sea con- 
tribution to it being quite essential. 

An important question is l/iVj, corrections to baryon observables. These 
can be classified in two groups: one comes from meson loops and is therefore 
accompanied by an additional small factor ^^ I/Stt^, the second arises from a 
more accurate account for the quantization of the zero rotational modes. The 
corrections of the second type are not accompanied by small loop factors, and 
may be quite substantial: after all in the real world A*"c = 3 so a 30% correction 
is not so small. Such corrections for certain quantities have been fished out 
in Refs. (63,64) in the two-flavor case and in Ref.(65) for three flavors. These 
corrections work in a welcome direction: they lower the fraction of nucleon 
spin carried by quark spins and increase the flavor non-singlet axial constants. 

Recent applications of the model are to parton distributions in nucleon, 
including the so-called skewed distributions, and to the wave function of the 
nucleon on the light cone. These topics deserve special sections. 

8 Nucleon structure functions 

The distribution of quarks, antiquarks and gluons, as measured in deep in- 
elastic scattering of leptons, provides us probably with the largest portion of 
quantitative information about strong interactions. Until recently only the 
evolution of the structure functions from a high value of the momentum trans- 
fer Q^ to even higher values has been successfully compared with the data. 
This is the field of perturbative QCD, and its success has been, historically, 
essential in establishing the validity of the QCD itself. However, the initial 
conditions for this evolution, namely the leading-twist distributions at a rela- 
tively low normalization point, belong to the field of non-perturbative QCD. If 
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we want to understand the vast amount of data on unpolarized and polarized 
structure functions we have to go into non-perturbative physics. 

The Chiral Quark-Sohton Model presents a non-perturbative approach to 
the nucleons, and it is worthwhile looking into the parton distributions it pre- 
dicts. Contrary to several models of nucleons on the market today, it is a 
relativistic field-theoretic model. This circumstance becomes of crucial impor- 
tance when one deals with parton distributions. It is only with a relativistic 
field-theoretic model one can preserve general properties of parton distribu- 
tions such as 

• relativistic invariance, 

• positivity of parton distributions, 

• partonic sum rules which hold in full QCD. 

There are two seemingly different ways to define parton distributions. The 
first, which we would call the Fritsch-Gell-Mann definition, is a nucleon matrix 
element of quark bilinears with a light-cone separation between the quark ip 
and ip operators. According to the second, which we would call the Feynman- 
Bjorken definition, parton distributions are given by the number of partons 
carrying a fraction x (the Bjorken variable) of the nuckxin momentum in the 
nucleon infinite-momentum frame. See Feynman's bookE£l for the discussion of 
both definitions. In perturbative QCD only the Fritsch-Gell-Mann definition 
has been exploited as it is very difficult to write down the nucleon wave function 
in the infinite-momentum frame, which is necessary for the Feynman-Bjorken 
definition. 

Despite the apparent difference in wording, it has been shown for the first 
time, within the field-theoretic Chiral Quark-Soliton Model, that the two defi- 
nitions are, in fact, equivalent and lead to identical working formulae for com- 
puting parton distributions: in Ref.(67) the first definition has been adopted 
while in Ref.(68) the second was used. The deep reason for that equivalence is 
that the main hypothesis of the Feynman-Bjorken parton model is satisfied in 
the Chiral Quark-Soliton Model, namely that parton transverse momenta do 
not grow with Q^. 

Let us point out some key findings of Refs. (67,68). 

8. 1 Classification of quark distributions in N^, 

Since the nucleon mass is 0{Nc) all parton distributions are actually functions 
of xNc- Combining this fact with the known large- TVc behavior of the integrals 
of the distributions over x one infers that all distributions can be divided into 
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'large' and 'small'. The 'large' distributions are, for example, the unpolarized 
singlet and polarized isovector distributions, which arc of the form 

i?i-^g°(x) ~ N^fixN,), (98) 

where /(y) is a stable function in the large- iVc limit. On the contrary, the 
polarized singlet and unpolarized isovector distributions give an example of 
'small' distributions, having the form 

D'^-'\x) ^ NJixN,). (99) 

One, indeed, observes in experiment that 'large' distributions are substantially 
larger than the 'small' ones. 

It should be mentioned that in the large- A^c limit there exist certain strict 
inequalities between various polarized and unpolarized distributions .Lj 

8.2 Antiquark distributions 

In the academic limit of a very weak mean pion field in the nucleon the Dirac 
continuum reduces to the free one (and should be subtracted to zero) while the 
bound-state level joins the upper Dirac continuum. In such a limit there are 
no antiquarks, while the distribution of quarks becomes q{x) — N^d{xNc — 1). 
In reality there is a non-trivial mean pion field which (j) creates a bound-state 
level, (ii) distorts the negative-energy Dirac continuum. As the result, the 
above (5-function is smeared significantly, and a non-zero antiquark distribution 
appears. 

An inevitable consequence of the relativistic invariance is that the bound- 
state level produces a negative contribution to the antiquark distribution. | 
The antiquark distribution becomes positive only when one includes the con- 
tribution of the Dirac continuum. Numerically, the antiquark distribution 
appears to be sizeable even at a low normalization point, in accordance with 
phenomenology. 

8.3 Sum rules 

The general sum rules holding iiiJiill QCD are automatically satisfied in the 
Chiral Quark-Soliton Model: inOHthe validity of the baryon number, isospin, 

^This is also true for any nucleon model with valence quarks, for example for any variant of 
bag models. Bag models are essentially non-relativistic, so they fail to resolve this paradox. 
In order to cure it, one has to take into account contributions to parton distributions from all 
degrees of freedom involved in binding the quarks in the nucleon. That can be consistently 
done only in a relativistic field-theoretic model, like the one under consideration. 
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total momentum and Bjorken sum rules has been checked. In fact, it is for 
the first time that nucleon parton distributions at a low normalization point 
have been consistently calculated in a relativistic model preserving all general 
properties. 

8.4 Smallness of the gluon distribution 

As many times stressed in this paper, the whole approach of the Chiral Quark- 
Soliton Model is based on the smallness of the algebraic parameter {MpY 
where 1/ p is the UV cutoff off the low-momenta theory, for example, the aver- 
age size of instantons in the vacuum. This p is the size of the constituent quark, 
while the size of the nucleon is, paramctrically, \/M . Computing parton distri- 
butions in the model one is restricted to momenta k ^ M k, 350 MeV, so that 
the internal structure of the constituent quarks remains unresolved. There are 
no gluons in the nucleon at this resolution scale; indeed, the momentum sum 
rule is satisfied with quarks and antiquarks only. 

However, when one moves to the resolution scale of 600 MeV or higher, 
the constituent quarks cease to be point-like, and that is at this scale that 
a non-zero gluon distribution emerges. Having a microscopic theory of how 
quarks get their dynamical masses one can, in principle, compute the non- 
perturbative gluon distribution in the constituent quarks. What can be said 
on general grounds is that the fraction of momentum carried by gluons is of 
the order of {MpY ~ 1/3, which seems to be the correct portion of gluons at 
a low normalization point of about 600 MeV where the normal perturbative 
evolution sets in. 

8.5 Comparison with phenomenology 

There are several parametrizations of the nucleon parton distributions at a 
relatively low normalization point, which, after their perturbative evolution to 
higher momentum transfer Q'^, fit well the numerous data on deep inelastic 
scattering. The paLametrization most convenient for our purpose is that of 
Gliick, Reya et al. l30 who pushed the normalization point for their distribu- 
tions to as low as 600 MeV, starting from the perturbative side. In Refs. (67,68) 
parton distributions following from the Chiral Quark-Soliton Model have been 
compared with those of Refs. (70,71), see Fig. 5. There seems to be a good 
qualitative agreement though the constituent quark and,^tiquark distribu- 
tions appear to be systematically 'harder' than those of.E^IO This deviation is 
to be expected since the calculations refer to even lower normalization point 
where the structure of the constituent quarks themselves is not been taken into 
account yet, see above. 
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Figure 5: Singlet unpolarized structure functions x{[u{x) + d{x)] ± [u{x) + d{x)]}/2 from 
Ref. (68) compared to the phenomenological parametrization of Ref. (70). The dashed Une 
corresponds to the sum and the soUd hne corresponds to the difference of quark and antiquark 
distributions. 



8. 6 Flavor asymmetry of the sea 

An important feature of the model is the large flavor asymmetry of the sea. 
Namely, isovector aniigwarfc distributions are not small, especially the polarized 
Au — Ad distribution which is of the leading order in Nc, while the unpolar- 
ized u — d distribution is of the subleading order. Phenomenologically, it is a 
welcome feature but which is hard to accommodate in other nucleon models. 
In the Chiral Quark-Soliton Model it arises in a most natural way since the 
nucleon is bound by the isospin 1 pion field. 

In Fig. 6 we present th&,calculated u — d distribution O as compared to 
the experimental data. SeeO for the latest discussion of the consequences of 
large flavor asymmetry of the sea on the description of the data. 

Finally, we mention that the Chiral Quark-Soliton Model has been used 
to make predictions on the transversity spin distributionso 



9 Skewed parton distributions 

Recently a new type of parton distributionsE3£J has attracted considerable in- 
terest, the so-called skewed (also called non-forward, off-forward, non-diagonal) 
parton distributions (the SPD's) which are generalizations of (i) usual parton 
distributions, {it) distribution amplitudes and {iii) of the elastic nucleon form 
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Figure 6: The antiquark asymmetry d{x) — u{x) in the proton computed in Ref. (72) and 
perturbatively evolved to Q = 7.35 GeV, as compared to the FNAL E866 data of Ref.(75). 



factors, in the case of vacuum quantum numbers in the t channel. At the same 
time there is a wide range of processes where one probes skewed densities which 
do not have a diagonal analog.E3 

The SPD's arc not accessible in standard inclusive measurements. How- 
ever, they can be measured in diffractive photQpr.nduction of the Z boson,L3 in 
diffractive electammduction of vector mesons jI3SEil in deeply- virtual Comp- 
ton scattering OoH and in phacd exclusive electroproduction of mesons at 
moderate values of Bjorken x.L2rE3 

The number of papers devoted to skewed distributions is rather large by 
now and continue to grow rapidly. Computation of the SPD's at a low normal- 
ization point in a realistic model of nucleons is of great importance. While the 
usual parton distributions have been measured in a variety of different experi- 
ments and can be confronted a posteriori with model calculations, in the case 
of the off-forward distributions the situation is opposite: model calculations of 
the distributions at a low normalization point are required to determine the 
very feasibility of measuring the SPD's. 

First roDdel calculations of skewed distributions were performed in the 
bag modcLEj However, bag models encounter severe problems in applications 
to parton distributions. It is usually assumed that the three quarks in the bag 
give rise only to quark distributions; however, they produce also an antiquark 
distribution with a negative sign, as an inevitable consequence of relativistic 
invariance. To overcome this problem one needs to add the contribution from 
forces which bind the quarks - in the case of the bag model it is the mysterious 
bag surface whose 'structure function' will hardly be ever defined. The Chiral 
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Quark-Soliton Model has no difficiiltjes of that kind, and skewed distributions 
have been calculated in the modeLa They appear to be qualitatively different 
from those of the bag model. We emphasize that this calculation is done in a 
self-consistent, relativistic-invariant model which ensures the positivity of both 
quark and antiquark distributions. 

Skewed parton distributions are defined through non-diagonal matrix ele- 



ments of the prodiii 
the notations of JiL 



3f quark fields at a light-cone separation. We shall use 
3 for the two functions depending on three variables: 



^e*^"(P'|V;(-An/2)V'(An/2)|P> =iJ(x,^,A2) {u{P')^u{P)) 

ZTT 

+ ^ ^(^' ^' ^') (^(^') *^'"'"A<^- "(^)) • (100) 

ZAnN 

Here n^ is a light-cone vector, n^ = 0, n ■ {P + P') = 2, and A is the four- 
momentum transfer, A = P' — P. A4n as usually denotes the nucleon mass, 
and u{P), u{P') are the Dirac spinors of the nucleon in the initial and final 
states. The off-forward quark distributions, H{x,£^,A'^) and E{x,^,A'^), are 
regarded as functions of the variable x, the square of the four-momentum 
transfer A^ and of its longitudinal component ^ = — (n • A) 

In the forward c ase, P = P', both A and ^ are zero, and the second term 
on the r.h.s. of Eq. (|lO0| ) vanishes. The function H becomes the usual parton 
distribution function, 

H{x,^^0,A^ ^0)^q{x). (101) 



On the other hand, taking the first moment of Eq. (|100|) one reduces the 
operator on the l.h.s. to the local vector current. The dependence of H and E 
on ^ disappears, and the functions reduce to the usual electric and magnetic 
form factors of the nucleon, 



J dxH{x,t^^) = ^i(A'), 
J dxE{x,^,A^) = F2(A2). 



The large Nc limit which is the basis of all calculations in the Chiral 
Quark-Soliton Model allows to calculate off-forward distributions in the region 
a;,^ ~ 1/Nc, A^ ~ N^. Also, it introduces a classification of the skewed 
distribution in Nc- in the leading order in Nc only isosinglet part of H{x, ^, A^) 
and isovector part of E{x,^,A'^) appear. Isovector H{x,£,,A'^) and isosinglet 
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E{x,£^,A'^) are subleading in Nc'. they appear after taking into accomit the 
finite angular velocity of the soliton rotation. 

The calculation of the off-forward parton distributions proceeds in much 
the same way as that of the usual parton distributions. A typical off-forward 
distribution obtained from the Chiral Quark-Soliton Model is plotted in Fig. 
7. Notice the oscillating character of the distribution with sham, crossovers 
near the points x — ±^/2. It can be shown on general grounds Eil that they 
should be present almost in any self-consistent field-theoretic model, n At small 
A (that is in the forward limit) these points divide the interval — 1 < x < 1 
into three regions: (i) x > ^/2, where the skewed distribution H{x,£,, A^) is 
close to the quark parton distribution, (ii) x < — C/2 where it is close to the 
antiquark parton distribution (with a minus sign) and {Hi) the intermediate 
region where it has a sharp crossover. The crossover is due to the fact that both 
quark and antiquark distributions should be positive. It can be shownEJ that 
the dependence of the quark mass M(k) on the quark virtuality converts a jump 
(with infinite derivative) in the skewed distributions into a sharp crossover with 
the width of the order '-^ (1/Mp)^. 

The fast crossover of H{x, ^, A^) at a; = ±^/2 may have interesting physical 
implications. For example, it may lead to a considerable enhancement of the 
deeply-virtual Compton scattering cross section: parametrically it is enhanced 
by a factor (logAfp)^. 

Applications of the Chiral Quark-Soliton Model to polarized skewed distri- 
butions can be found in Ref.(85). The model predicts a_cQiisiderable enhance- 
ment of the hard exclusive production of charged pionsEj'cIl accompanied by a 
large azimuthal spin asymmetry.E3 Also, it appears possible to obtain certain 
model-independent relations between different skewed distributions (based on 
the large Nc limit), to calculate the transition iV — > A distributions, etc. 

10 Light-cone wave function of the nucleon 

Probably, the most complete information about processes at high energies in- 
volving nucleons can be extracted from its wave function on the light cone or, 
else, in the infinite momentum frame. Using the infinite momentum frame al- 
lows one to separate the wave function of the nucleon from vacuum fluctuations. 
The wave function is, generally, a Fock vector which determines amplitudes to 
find inside the nucleon three quarks, or three quarks plus a quark-antiquark 
pair, or three quarks plus one gluon, etc. The wave function squared gives the 
probability to find quarks inside the nucleon with given fractions of nucleon 
longitudinal momentum and some transverse momenta. Integrating it over 

''They are absent in bag model calculations £j 
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Figure 7: Prediction for the isosinglet skewed distribution H(x,^, A'^) for A|, = (A|, = 
— A^ — g^M|^) and g = 0.3. Dashed line: contribution from the discrete level; dashed-dotted 
line: contribution from the Dirac continuum; solid line: the total distribution (sum of the 
two). The vertical lines mark the crossover points x = ±5/2. 



transverse momenta and summing up over all possible configurations with one 
quark (or antiquark) fixed one arrives to standard parton distributions. 

By the wave function in a 'narrow sense' (j){zi, Z2, z^) (else called distribu- 
tion amplitude) people usually understand the amplitude to find exactly three 
quarks inside the nucleon with given fractions zi,Z2,z^ of the nucleon longi- 
tudinal momentum. It determines the high-energy asymptotics of exclusive 
processes such as the asymptotics of electromagnetic formfactors, tke decays 
J ftp -^ NN, etc. It is rather well known from experiment (see, e.g. ESl). 

The Chiral Quark-Soliton Model presents a unique possibility to determine 
the whole Fock state vector of the nucleon on the light cone. The nucleon 
wave function (as well as parton distributions) are subject to the evolution 
in perturbation theory, therefore we are actually speaking about the wave 
function at a low normalization point. This work is in progress, so we present 
some preliminary results.E^I 

In the Chiral Quark-Soliton Model it is easy to define the wave function 
of the nucleon in the rest frame. Indeed, the nucleon in the model represents 
quarks in the Hartree approximation in the self-consistent pion field. In this 
approximation the full wave function is the product of one-particle states. 



Yl /d3x/(")(a;)*+(x)|0>. 



\^N>= 11 ld'xf(^>{x)^+{x)\0>, (102) 

occupied ^ 

where f"{x) are the eigenf unctions of the time-independent Dirac operator in 
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the external pion field, Hf^"'> = £■(")/("), with the hamiltonian @. The 
product goes over all occupied levels, i.e. over the discrete 'valence' level and 
over the negative-energy continuum. 

Actually we need the wave function not in terms of the Dirac continuum 
but rather in terms of quarks and antiquarks. It is well known that the wave 



function (102) corresponds to the so-called coherent exponent: 



color ' 



xexp f dkidk2a^' + {ki)ex,xAki,k2)h^'+ik2), (103) 

where a+,b+ are quark and antiquark creation operators. It represents Nc 
quarks with the wave function F^{k) and a number of additional quark- 
antiquark pairs whose wave functions in the external pion field are 0AiA2 (^ii fe); 
they can be calculated using the Feynman Green function at finite time, which 
can be constructed from the Dirac equation. As to the single-quark wave func- 
tion F^{k) it gets contributions both from the discrete level and from the Dirac 
continuum described by the pair wave function Qx^x^iki, fc2). 
We next proceed as follows: 



The wave function (103) is translated to the infinite momentum frame. 



This is possible as the model is relativistically invariant. 

• The pair wave functions are calculated using the technique of interpo- 
lating formulae described in subsection 4.1. It is interesting to note that 
higher terms of this expansion. give the wave functions with 1,2,. . . addi- 
tional pions on the light cone.E3E2l This is to be expected, as pions are 
the only agents inducing the interaction in the model. 

• The result should be averaged with the rotational wave function defined 
in section 6.2 to project the general wave function onto a given rotational 
state of the soliton - the nucleon, the A, etc. 

The three-quark contribution to the (proton, spin up) wave function is: 



|$jv, proton T>= $(2:1, pix,Z2,P2±,Z3,P3±) WUl) > \ui{2) > |dT(3) > . 

(104) 
After integrating over transverse momenta of quarks one gets the distribution 
amplitude. The discrete-level contribution to it has the form: 
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Figure 8: Nucleon wave function </>(zi, Z2, 1 — 21 — 22)- Only the 'valence' level contribution 
is included. 



0(2^1,22,23)= / (fx_L{g{zi,xi_)g{z2,x±_)g{z3,x±) 
+ x\j{z2,xi_) [25(zi,a;_L)j(23,a:^±) - i{zi,Xi_)g{z^,XA_)]] , (105) 



where 

g{z,XA_) = - 
i{z,xi_) 



3/2 



2 

^3/2 



d.T3 CXp(-i(zAlAr - i^lcvcl)) 

dx'i exp(-i(zA^7v - -Eicvoi)) 



h{r) j{r) 

r 



-j(0 



,(106) 
(107) 



Here h{r) and j{r) are the wave functions of the discrete level (p3) and ii^icvci 
is its energy. 



Let us note that the second term in curly brackets in Eq. (105) introduces 



an asymmetry of the wave function under the exchange of z^'s. This effect is 
due to the finiteness of Nc and appears as a result of averaging over nucleon 
rotational wave function. At Nc — > 00 the wave function of the nucleon would 
be completely symmetric. 



The calculated distribution amplitude (105) is presented in Fig. 8. Though 
not fully symmetric under the interchange of Zi^s, it is much closer to the 
asymptotic nucleon wave function, 



6"^y'°P*(zi, Z2, 23) = 120 Z1Z2Z3 6{ZI +Z2+Z3~ 1), 
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(108) 



than the wave function suggested by Chernyak, Oglobhn and ZhitniskyE3 many 
years ago on the basis of. their analysis of the QCD sum rules. Fl It seems to be 
supported by the data.E3 



11 Conclusions 

The Chiral Quark-Soliton Model of baryons is a simple but still a non-trivial 
reduction of the full-scale QCD at low energies. It emphasizes the spontaneous 
chiral symmetry breaking in QCD, accompanied by the appearance of the 
dynamical (or constituent) quark mass. We prefer the word 'dynamical': first, 
because it is, indeed, dynamically generated, second, because it is momentum- 
dependent. 

The momentum dependence of the dynamical quark mass M{k) is the 
key to understanding why the notion of constituent quarks have worked so 
remarkably well over 30 years in hadron physics. The point is, the scale A at 
which the function M{k) falls off appears to be much larger than M(0); the 
former determines the size of constituent quarks while the latter determines 
the size of hadrons. As a matter of fact these two distinctive scales come neatly 
from instantons. 

The Chiral Quark-Soliton Model fully exploits the existence of the two 
distinctive scales: it is because of them it makes sense to restrict oneself to 
just two degrees of freedom in the nucleon problem, namely, to massless or 
nearly massless (pseudo) Goldstone pions and to the constituent quarks with a 
momentum-dependent dynamical mass M{k). The scale A actually plays the 
role of the physical ultra-violet cutoff for the low-energy theory; its domain 
of applicability is thus limited to the range of momenta fc ~ M < A. This is 
precisely the domain of interest for the nucleon binding problem. 

A technical tool simplifying considerably the nucleon problem is the use of 
the large Nc logic. At large Nc the nucleon is heavy, and one can speak of the 
classical self-consistent pion field binding the Nc valence quarks of the nucleon 
together. The classical pion field (the soliton) is found from minimizing the 
energy of the bound-state level plus the aggregate energy of the lower Dirac 
continuum in a trial pion field. The 'valence' quarks sitting on the bound-state 
level appear to be strongly bound by the classical pion field. The system is 
relativistic. 

By quantizing the slow rotations of the soliton field in flavor and ordinary 
spaces one gets baryon states which are rotational excitations of the static 
'classical nucleon'. The classification of the rotational excitations depends on 



'The pirn wave function at a low normalization point as calculated from the instanton 
vacuum Eij also appears to be rather close to its asymptotic form. 
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the symmetry properties of the soUton field, but not on the details of dynamics. 
Taking the hedgehog ansatz one gets the following lowest baryon multiplets: 
octet with spin 1/2, decuplet with spin 3/2 (these are, indeed, the lowest mul- 
tiplets observed in nature) and antidecuplet with spin 1/2. This last multiplet 
contain baryons with exotic quantum numbers (in the sense that they cannot 
be composed of only three quarks); some of them are predicted to be relatively 
light and narrow resonances, and it would be of great interest to search for 
such states. 

By saying that all lightest baryons are nothing but rotational excitations 
of the same object, the 'classical nucleon', one gets many relations between 
members of baryon multiplets, which are realized with astonishing accuracy in 
nature. Especially successful are predictions which do not depend on dynami- 
cal quantities (like the values of moments of inertia) but follow from symme- 
try considerations only, and are therefore shared, e.g., by the Skyrme model. 
Predictions of the Chiral Quark-Soliton Model which do depend on concrete 
dynamics are, in general, also in good accordance with reality: the typical 
accuracy for numerous baryon observables computed in the model is about 
15-20%, coinciding with the expected theoretical accuracy of the model. To 
get a better accuracy one needs a better understanding of the underlying QCD 
vacuum and of the resulting effective low-energy theory. 

To our knowledge, the Chiral Quark-Soliton Model is today the only rel- 
ativistic field-theoretic model of the nucleon, and this advantage of the model 
becomes crucial when one turns from static characteristics of the nucleon to the 
numerous parton distributions. It is impossible to get a consistent description 
of parton distributions satisfying positivity and sum rules restrictions, with- 
out having a relativistic theory at hand and without taking into account the 
complete set of forces which bind quarks together. The leading-twist parton 
distributions calculated in the Chiral Quark-Soliton Model refer to a very low 
normalization point where the structure of the constituent quarks is not re- 
solved yet. They seem to be in qualitative agreement with paramctrizations 
of the DIS data at low Q^ though, not unnaturally, they appear to be more 
'hard'. One of the striking predictions of the model is large flavor asymmetry 
of antiquark distributions, which is favored by recent experiments. An even 
larger asymmetry is predicted for polarized antiquark distributions, but it has 
not been directly measured so far. 

The Chiral Quark-Soliton Model has been recently used to calculate some 
of the off-diagonal or skewed distribution functions which provide even a more 
detailed information on the nucleon structure than the usual (diagonal) parton 
distributions. It predicts a drastic oscillating form of the skewed distribution. 
Preliminary results have been obtained for the nucleon wave function on the 
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light cone. 

An impressive amount of data on parton distributions has been collected 
in the past years, and its amount will grow with new types of distributions 
becoming experimentally accessible. We know how to (perturbatively) evolve 
parton distributions from high to still higher values of Q^ but we do not under- 
stand how to explain the initial conditions for that evolution, since it involves 
truly nonperturbative physics. This is a real challenge to the theory of strong 
interactions. At present the only nonperturbative relativistic field-theoretic 
model of QCD at low energies is the 'Nucleons as Chiral Solitons' model. 
Therefore, we think that it will play a major role in understanding, describing 
and predicting data in the years to come. 
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